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ABSTRACT 



We present new radial velocities for 289 globular clusters around NGC 4636, the southernmost giant elliptical galaxy of the Virgo 
cluster. The data were obtained with FORS2/MXU at the Very Large Telescope. Together with data analysed in an earlier study 
(Schuberth et al. 2006), we now have a sample of 460 globular cluster velocities out to a radius of 12 arcmin (60 kpc) available - 
one of the largest of its kind. This new data set also provides a much more complete angular coverage. Moreover, we present new 
kinematical data of the inner stellar population of NGC 4636. We perform an updated Jeans analysis, using both stellar and GC data, 
to better constrain the dark halo properties. We find a stellar M/L-ratio of 5.8 in the R-band, higher than expected from single stellar 
population synthesis. We model the dark halo by cored and cuspy analytical halo profiles and consider different anisotropics for the 
tracer populations. Properties of NFW halos lie well within the expected range of cosmological simulations. Cored halos give central 
dark matter densities, which are typical for elliptical galaxies of NGC 4636's luminosity. The surface densities of the dark matter 
halos are higher than those of spiral galaxies. We compare the predictions of Modified Newtonian Dynamics with the derived halo 
properties and find satisfactory agreement. NGC 4636 therefore falls onto the baryonic Tully-Fisher relation for spiral galaxies. The 
comparison with the X-ray mass profile of Johnson et al. (2009) reveals satisfactory agreement only if the abundance gradient of hot 
plasma has been taken into account. This might indicate a general bias towards higher masses for X-ray based mass profiles in all 
systems, including galaxy clusters, with strong abundance gradients. 

Key words, galaxies: elliptical and lenticular, cD — galaxies: kinematics and dynamics — galaxies:individual: NGC 4636 



1. Introduction 

NGC 4636 is a remarkable elliptical galaxy. Situated at the 
Southern border of the Virgo galaxy cluster, and thus not in a 
very dense environment, its globular cluster system (GCS) ex- 
hibits a richness, which one does not find in othe r galaxies of 
comparable luminosity in a similar environments (iKissler et al.l 
Il994t lDirschetal.l[2005l) . Therefore, NGC 4636 offers the op- 
portunity to employ globular clusters (GCs) as dynamical trac- 
ers to investigate its dar k halo out to large r adii, which is rarely 
given. In our first study (ISchuberth et al.l2006l hereafter Paper I), 
we measured 174 GC radial velocities to confirm the exis- 
tence of a dark halo, previously indi cated by X-ray analyses 
(e.g. lLoewenstein & Mu shotzkv 2003) and tried to constrain its 
mass. Paper I also gives a summary of the numerous works re- 
lated to NGC 4636 until 2006, which we do not want to repeat 
here. iDirsch et al.l (|2005) presented a wide-field photometry in 
the Washington system of the globular cluster system of NGC 
4636, which was the photometric base for Paper I as well as 
for the present work. Also there, the interested reader will find 
a summary of earlier works. Noteworthy peculiarities of NGC 
4636 include the appearance of the supernova 1937 A (a bona- 



fide la event, indicative of presenc e of an intermedia t e-age pop- 
ulation), the high FIR-emission dTemi et alJ I2003L 120071), the 



chaotic X-ray fe atures in the inner region djones et al.l l2002t 
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Baldi et al.l2009h . pointing to feed-back effects from supernovae 
or enhanced nuclear activity in the past. TableQ] summarises the 
basic parameters related to NGC 4636. 

Since NGC 4636 is X-ray bright, it has be en the target of 
nume r ous X-ray studies (for earlier work s ee Matsushita et al.l 
11998b iJones et al.l l2002t Loewenstein & Mushotzkvl 120031 and 
references therein), offering the possibility to compare X-ray 
based mass profiles with stellar dynamical mass profiles, using a 
data base, which is not found elsewhere for a non-central ellipti- 
cal galaxy. 

The most recen t X-ray based analysis of the mass profile is 
from Joh nson et alJ (120091) whose work is based on Chandra X- 
ray data. The new feature in their work is the inclusion of an 
abundance gradient in the X-ray gas. We will show that, if this 
gradient is accounted for, the resulting X-ray mass profile is 
in good agreement with the one derived from our GC analysis, 
which may indicate the need to respect abundance gradients, if 
present, in any X-ray analysis. 

Two recent p ublica tions on th e NGC 4636 GCs are from 
iPark et all (120101) and iLee et al] (120101) . who present about 
100 new Subaru spectra of NGC 4636 GCs a nd discuss their 
kinematics in combination with the data from Schuberth et al.l 
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(2006). For the sake of brevity, we postpone a detailed compari- 
son of the different samples to a later publication. 

Our data from Paper I have also been used by 
IChakrabartv & Ravc haudhurv (2008) who employed their 
own code to obtain a mass profile. 

In the analysis presented in Paper I, the most important 
source of uncertainty is the sparse data at large radii: The dis- 
persion value derived for the outermost radial bin changes dras- 
tically depending on whether two GCs with extreme velocities 
are discarded or not. If these data are included, the estimate of 
the mass enclosed within 30 kpc goes up by a factor of ~ 1 .4 and 
the inferred dark halo has an extremely large (^ 100 kpc) scale 
radius. 

Moreover, the data presented in Paper I have a very patchy 
angular coverage since the observed fields were predominantly 
placed along the photometric major axis of NGC 4636 (see 
Fig. |U right panel). To achieve a more complete angular cov- 
erage and to better constrain the enclosed mass and the shape of 
the NGC 4636 dark halo, we obtained more VLT FORS 2 MXU 
data, using the same instrumental setup as in our previous study. 

In our present study, an important difference to Paper I is the 
distance, which is of paramount importance in the dynamical 
discussion. In Paper I, we adopted a distance of 15 Mp c, based on 
surfac e brightness fluctuation (SBF) measurements dTonrv et~ai] 
1200 lb . but already remarked that we consider this value to be a 
lower limit. A more recent re-cal ibration of the S BF brings the 
galaxy even closer to 13.6 Mpc (I Jensen et al.ll2003l) . 

On the other hand, the method of globul ar luminosity func- 
tions revealed a distance of 17.5 Mpc (Kis sler et alj 1 19941 
iDirsch et al.ll2005| . We do not claim a superiority of GCLFs 
over SBFs, but one cannot ignore the odd findings which re- 
sult from adopting the short distance. The specific frequency of 
GC s would assume a v alue rivalling that of central cluster galax- 
ies ([pirsch et aEl l2005l) . Moreover, the stellar M/L for which in 
PapeiT we already adopted a very high value of 6.8 in the R- 
band, would climb to 7.5 for a distance of 13.6 Mpc, which al- 
most dou bles the value expecte d for an old metal-rich popula- 
tion (e.g. Cap pellari et al. 2006). This is under the assumption 
that the total mass in the central region, for which we present 
independent kinematical data, is dominated by the stellar mass 
(the 'maximal disk' assumption). The alternative, namely that 
NGC 4636 is dark matter dominated even in its centre, would be 
intriguing, but would make NGC 4636 a unique case. To adopt a 
larger distance is the cheapest solution. The explanation for the 
discrepancy between SBF and GCLF distance might lie in some 
effect enhancing the fluctuation signal, either by an intermediate- 
age population (which in turn would require a lower M/L) or an- 
other small-scale structure, which went so far undetected. In the 
following, we adopt a distance of 17.5 Mpc. Thus, l'is ~5.1 kpc 
andl"=85pc. 

This paper is organised as follows: In Sect. [2] we give a brief 
description of the spectroscopic observations and the data re- 
duction. In Sect. [3] we present the combined data set, before dis- 
cussing the spatial distribution and the photometric properties of 
the GC sample in Sect. [4] In Section|5j we define the subsamples 
used for the dynamical study. The line-of-sight velocity distri- 
butions of the different subsamples are presented in Sect. |6] In 
Sect. |7] the subsamples are tested for rotation, and in Sect. [8] 
the velocity dispersion profiles are shown. In Sect. |9]we present 
new stellar kinematics for NGC 4636 and in Sect. [10] we sum- 
marise the theoretical framework for the Jeans modelling. The 
mass profiles are shown in Sect. 11. Sections 12 and 13 give the 
discussion and the conclusions, respectively. 



Table 1. NGC 4636 basic data 



NGC 4636 basic data 



Other names 


UGC 07878, VCC 1939 


NED 


Position (2000) 


12M2?49 S .8 +02°41'16" 


NED 




190.707792 2.687778 




Galactic coordinates 


/ = 297.75° b = 65.47° 




Galactic extinction 


A B = 0.050 


(1) 




A B = 0.118 


(2) 


Distance 


D = 17.5 Mpc 


(Sect. 1) 


Scale 


l"=85pc 1'= 5.09 kpc 




eta tipp moHnliiQ 


( m _ M\ — 31 ?? 

yffi lYi j — ji.z.^ 




Heliocentric velocity 


tViio = 906 ± 7kms~ l 


(3) 


Hubble type 


E0+ 


R3C 


Ellipticity 


e = 0.15 


(4) 


Position angle 


PA = 150° 


R3C 


dllCCLlVC latilUS 


d _ i m"7 (— 8 f\A \cr\r \ 


w 


Age 


1 J.J ± J.O v jyl 


PJ 


Metallicity 


7 — n mi -l n nnA 

Z, — U.UZJ ± U.UUD 


{.->) 


a /Fe 


[Li/rc] — \j.jLy it u.uu 


(J) 


Total blue mag 


b T = 9.78 


(6) 


Absolute blue mag 


B T = -21.05 




Colours 


U-B = 0.50 


(6) 




B-y = 0.87 






V-R = 0.67 






V-l = 1.30 




Stellar pop M/L 


M/L, = 3.74 


(7) 


Dynamical M/L 


M/L B = 12.2 


(8) 


X-ray luminosity 


L x = 1.78±0.10xl0 41 ergs/s 


(9) 


Nuclear X-ray emission 


< 2.7 x 10 38 ergs/s 


(10) 


Central black hole 


M S mbh~8x 10 7 M o 


(11) 



Note 1. UGC = Uppsala General Catalogue of Galaxies 
iNilsorJ 119731) . VCC = Virgo Cluster Catalogue (Bingee li etaD 
119851), NED = N ASA/IPAC Extragalactic Database 
(http : //nedwww. ipac . caltech . ed u), R3C = Third R e ferenc e 
Catalogue of Bright Galaxies (Ide Vaucoule urs et al J 1 1 99 ll) 
References: (D lBurstein & Heilesl dl982h: (2) ISchleeel et"al 



(1998j); (3) Paper I; (4) iBenderet 
d2007T> : ( 6) iPrugniel & Heraudeauf 
2003); (8) iKronawitter et all feOOO) 



( 



.19941) : (5) lAnnibali et al 
19981): (7) dMaraston et al. 
(9) iForman et all dl985l) : (1Q) 



Loewenstein et al.ld200lh : (ll) lMerritt & Ferraresel(l200lh 



2. Observations and data reduction 

The data were acquired using the same instrumental setup as 
described in Paper I. Therefore, we just give a brief description 
here. 



2.1. VLT-observations 

The observations were carried out in service mode at the 
European Southern Observatory (ESO) Very Large Telescope 
(VLT) facility on Cerro Paranal, Chile. We used the FORS 2 
(FOcal Reducer/low dispersion Spectrograph) equipped with 
the Mask EXchance Unit (MXU). The programme ID is 
075.B-0762(B). 

Pre-imaging of twelve fields, shown in the left panel of 
Fig.Q] was obtained in April 2005. The GC candidates were se - 
lected from the photometric catalogue by IDirsch et al.l (I2005I) . 
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Fig. 1. NGC4636 GC spectroscopic data. The 33' x 33' DSS image shown in both panels is centred on NGC4636, North is to the 
top and East is to the left. Left: Location of the new fields. Note that only the inner 4 .'3 x 6.'8 of the FORS 2 field-of-view are 
shown, since the positions of the slits are confined to this area. Right: All NGC4636 velocity-confirmed GCs. Squares show the 
data set from Paper I, circles the new data. The white dashed circles have radii of 2.'0, 8.'0 and 13.'5, corresponding to approximately 
10, 41, and 69 kpc (at D = 17.5 Mpc), respectively. 



The targets have colours in the range 0.9 < C — R < 2.1 and 
R-bmd magnitudes brighter than 22.5 mag. 

The spectroscopic masks were designed using the ESO 
FORS Instrumental Mask Simulator (FIMS) software. In con- 
trast to Paper I, where object and sky spectra were obtained from 
different slits, we chose to observe the GC candidates through 
longer slits and to extract the background from the same slit. The 
slits for point-sources have a width of 1 " and a length of 4". The 
spectroscopic observations were carried out in service mode in 
the period May 2 nd to June 7 th , 2005. We used the Grism600B 
which gives a resolution of ~ 3 A. A total of twelve masks were 
observed with exposure times between 3600 and 5400 seconds. 
To minimise the contamination by cosmic-ray hits, the observa- 
tion of each mask was divided into two or three exposures. A 
summary of the MXU observations is given in Table|2] 

2.2. Data reduction 

Prior to bias subtraction, the f smosaic-script, which is part of 
the FIMS-software was used to merge the two CCD-exposures 
of all science and calibration frames. The spectra were traced 
and extracted using the Iraf apall package. 

For the extraction of point-sources, we chose an aperture 
size of 3 (binned) pixels, corresponding to 0.75". In each slit, 
we interactively defined an emission-free background region, 
the pixel values of which were averaged and subtracted off the 
spectrum during the extraction. This sky-subtraction was found 
to work very well, only the strongest atmospheric emission lines 
left residuals. The spectra were traced using the interactive mode 
of apall, and the curve fit to the trace was a Chebyshev polyno- 
mial of order 3-11. The use of this wide range of polynomials is 
motivated by the fact that the characteristics of the bottom CCD 
('slave') deviate from those of top ('master') CCD. Besides hav- 



ing a worse point spread function, the traces of spectra on the 
'slave' CCD are more contorted. Thus, the tracing of spectra 
on the 'slave' CCD required the use of the higher-order poly- 
nomials. On the 'master' frame, an order of 3-5 proved to be 
sufficient. 

The wavelength calibration is based on Hg-Cd-He arc-lamp 
exposures obtained as part of the standard calibration plan. The 
one-dimensional arc spectra were calibrated using the Iraf task 
identify. Typically, 17-22 lines were identified per spectrum 
and the dispersion solution was approximated by a 5 th -order 
Chebyshev polynomial. The rms-residuals of these fits were 
about 0.05 A. 

The wavelength calibrations were obtained during day-time, 
after re-inserting the masks and with the telescope pointed at 
zenith. Consequently, the offsets introduced by instrument flex- 
ure and more importantly, the finite re-positioning accuracy of 
MXU have to be compensated for. To derive the corresponding 
velocity corrections we proceeded as follows: The wavelength 
of the strong [OI] 5577 A atmospheric emission line was mea- 
sured in all raw spectra of a given mask as a function of the 
slit's location on the CCD (orthogonal to the dispersion axis). 
The correction for each aperture was determined from a 2 nd or- 
der polynomial fit to the wavelength-position data. The average 
magnitude of this correction was ~ 30kms . On a given mask, 
the wavelength drift between the lower edge of the slave chip to 
the top edge of the master chip corresponded to velocity differ- 
ences between 10 and 90 km s . 

2.3. Velocity determination 

To determine the radial velocities of our spectroscopic tar- 
gets, we proceeded as described in Paper! We used the 
same spectrum of the dwarf elliptical galaxy NGC1396 
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Table 2. Summary of NGC 4636 VLT FORS2/MXU observations 



Mask ID 


Centre Position 
(J 2000) 


Obs. Date 


Seeing 


# Exp. 


Exp. Time 

(sec) 


# Slits 


r Ul 


1 9 -zn-04 


+09-36-90 


ZUUJ-ID-U4 


V. Id 




J4UU 


J / 


r uz 


12-43-03 


+02-40-46 


ZUUJ-UJ-UJ 


ft" on 


D 




O / 


r Uj 


12:43:04.0 


+02:47:07.0 


ZUUj-Uj-UZ 


t. uu 






OJ 


In C\A 

r U4 


12:42:47.0 


+02:49:18.0 


ZUUJ-UJ-U / 


u. yu 


z 




A 1 
04 


F05 


12:42:45.0 


+02:43:37.0 


2005-05-05 


a/ 7i 


3 


5400 


63 


F06 


12:42:54.0 


+02:35:28.0 


2005-05-31 


a' 71 


2 


3600 


71 


F07 


12:42:39.0 


+02:41:16.0 


2005-05-31 


a/ 90 


2 


3600 


75 


F08 


12:42:30.0 


+02:38:35.0 


2005-05-05 


a/ 66 


2 


3600 


67 


F09 


12:42:24.0 


+02:47:22.0 


2005-05-05 


a/ 63 


3 


5100 


63 


F10 


12:42:23.0 


+02:42:35.2 


2005-06-01 


a/58 


3 


5400 


68 


F12 


12:43:14.0 


+02:40:46.0 


2005-06-02 


a/ so 


3 


5400 


62 


F13 


12:42:34.0 


+02:49:07.0 


2005-06-07 


0?76 


3 


5400 


68 



Note 2. ESO program ID 75.B-0762(B). The seeing values were measured from the acquisition images obtained prior to the corresponding 
spectroscopic observations. 




Fig. 2. Velocity uncertainties. Left: Comparison of the GC velocity estimates obtained with the two different template spectra. 
For graphic convenience, the error-bars (the two /rcor-uncertainties added in quadrature) are only shown for objects where the 
estimates differ by more than 75kms~'. The grey area shows the region encompassing 68 per cent of the data points. Middle: 
Velocity uncertainties as determined by fxcor. Left sub-panel: fxcor-uncertainty vs. /^-magnitude. Crosses and dots represent 
blue and red GCs, respectively. The solid and dashed curve show a quadratic fit to the blue and red GC data, respectively. Right 
sub-panel: fxcor-uncertainty vs. colour. The dashed line at C—R = 1.55 indicates the colour dividing blue from red GCs. The 
solid line illustrating the increase of the uncertainties towards bluer colours is a 2nd order polynomial fit to the data. The dashed 
lines at Av = 65kms~ I show the limit adopted for the sample definition in Sect. 15.21 Right: Comparison of the duplicate velocity 
measurements. Dots and crosses represent GCs and stars, respectively. The dashed lines show the r.m.s. scatter of 84kms~ 1 found 
for the GCs. 



(v = 8 15 ± 8 km s '. lDirsch et al.ll2004l) as template, and the ve- 
locities were measured using the Iraf fxcor task, whic h is based 
upon the technique described in lTonrv & Davisl(ll979l) . 

To double check our results, we also used the spectrum of 
one of the brightest globular clusters found in this new dataset 
(object f 12-24, v helio = 980 + 15 kms" 1 , m R = 19.9 mag, C-R = 
1.62) as a template. Note that these are the sam e templates we 
used in our recent study of the NGC 1399 GCS ISchuberth et all 
d2009l) . 

The wavelength range used for the correlation was AA 4100- 
5180 A. The blue limit was chosen to be well within the domain 
of the wavelength calibration. Towards the red, we avoid the 
residuals from the relatively weak telluric nitrogen emission line 



at 5199 A. In the few spectra affected by cosmic ray residuals or 
bad pixels, the wavelength range was adjusted interactively. 

The left panel of Fig. [2] shows the difference of the veloc- 
ities determined using the the 'GC and the 'galaxy' template 
as function of /^-magnitude. As expected, the scatter increases 
for fainter magnitudes. The r.m.s of the residuals is 35kms~', 
which is comparable to the mean velocity uncertainties returned 

by fxcor (Av gc = 37, Av ga i = 48kms~ 1 ). For each spectrum, 
we adopt the velocity with the smaller fxcor-uncertainty. 

The middle panel of Fig. [2] shows the fxcor-uncertainties 
as function of /^-magnitude (left sub-panel). The solid (dashed) 
curve is a quadratic fit to the blue (red) GCs. While the uncer- 
tainties increase for fainter objects, one also notes that, at a given 
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magnitude, the blue GCs, on average, have larger velocity un- 
certainties. The same trend can be seen in the right sub-panel, 
where we plot colour-dependence of the uncertainties: The bluer 
the objects, the larger the uncertainties. This trend is likely due to 
the paucity of spectral features in the spectra of bluest, i.e. most 
metal-poor GCs. 

3. The combined data set 

In this section, we describe the database used for the dynamical 
analysis. The final catalogue as given in Table lD . 2l combines the 
new velocities with the data presented in Paper I and the pho- 
tometry by Dirsch et al. (200 51, D+05 hereafter). 

Our data base has 893 entries, (547 new spectra, 346 veloc- 
ities from the catalogues presented in PapeiT). The new spectra 
come from 463 unique objects (327 GCs and 136 Galactic fore- 
ground stars3). Of these, 289 GCs and 116 stars were not tar- 
geted in the previous study. The spectra of background galaxies 
were discarded from the catalogue at the stage of the fxcor ve- 
locity determination, and there were no ambiguous cases. 
Thus, including the 171 (170) unique GCs (foreground stars), 
from Paper I, the data set used in this work comprises the veloc- 
ities of 460 individual bona-fide NGC4636 GCs and 286 fore- 
ground stars. 

3. 1 . Photometry database 

As in Paper I, we used the Washington photometry by D+05 to 
assign C-R colours and /?-band magnitudes to the spectroscopic 
targets. 

For some objects with velocities in the range expected for 
GCs, however, no photometric counterpart was found in the fi- 
nal photometric catalogue published by D+05. The reason for 
this is that the D+05 catalogue lists the point sources in the field , 
which were selected using the DAOPhotll dStetsonlll987lll992l) 

and 'sharpness' parameters. While this selection, as desired, 
rejects extended background galaxies, it also removes those GCs 
whose images deviate from point sources. The GCs which were 
culled from the photometric data set fall into two categories, the 
first encompassing objects whose images have been distorted be- 
cause of their location near a gap of the CCD mosaic, a detector 
defect (e.g. a bad row) or a region affected by a saturated star. 
The second group are GCs whose actual sizes are large enough to 
lead to slightly extended images on the MOSAIC images which 
have a seeing of about l'.'O. 

Using the 'raw' version of the photometry database, we 
were able to recover colours and magnitudes for 154 objects, 
80 of which are GCs (according to the criteria defined below in 
Sect. 13.3b . As can be seen from the left and middle panels of 
Fig. [5] many of the brightest GCs were excluded from the D+05 
catalogue. 

3.2. Duplicate measurements 

In our catalogue, we identify 131 (94) duplicate, 8 (7) triple mea- 
surements, and 607 (359) objects were measured only once (the 
corresponding numbers for the GCs are given in parenthesis). 
For objects with multiple measurements we combine the mea- 
surements using the velocity uncertainties as weights. The ve- 
locities from the duplicate measurements are compared in the 
right panel of Fig. [2] The r.m.s. scatter of 84km s found for 

1 Following Paper I, we consider objects with v hc i io below 350 km s~ 1 
as foreground stars (see Sect. 13.3) 



the GCs is about twice the average velocity uncertainty quoted 
by the cross-correlation programme. Random offsets in velocity 
might be introduced by targets which are not centred in the slit 
in the direction of the dispersion direction. 

3.3. Separating GCs from Galactic foreground stars 

In the left panel of Fig. [3] we plot the C-R colour versus helio- 
centric velocity: The GCS of NGC 4636 occupies a well-defined 
area in the velocity-colour plane. The foreground stars, cluster- 
ing around Zero velocity, have a much broader colour distribu- 
tiotfl 

As in PapeiT, we adopt v'heiio = 350kms _I as lower limit 
for bona-fide NGC 4636 GCs. Although both classes of objects 
are well separated in velocity space, there is a number of uncer- 
tain cases: The seven objects with velocities in the range 300- 
400 km s (indicated by the grey areas in both panels of Fig. [3} 
merit closer scrutiny. These objects are discussed in Sect. 15.1.21 
where we remove the likely outliers from our sample. 

The following section gives a description of the photomet- 
ric properties and the spatial distribution of the 460 bona-fide 
NGC 4636 GCs in our velocity catalogue. 

4. Properties of the GC sample 

Our kinematic sample now comprises a total of 460 GCs, which 
is more than 2.5 times the number of GCs used in Paperl. After 
NGC 1399 dRichtler et al l [20041 l200l [Schuberth et al. 2009) 
with almost 700 GC velocities, Cen A dWoodlev et al) 12010) 
with about 560 velocities, and M87 (?) with 488 velocities, this 
is the fourth-largest GC velocity sample to date. 

4.1. Spatial distribution 

The spatial distribution of the GCs with velocity measurements 
is shown in Fig. |4] The left panel plots the distribution of the 
galactocentric distances of the GCs as unfilled histogram. The 
azimuthal distribution is shown in the right panel. Compared to 
Paperl (grey histogram), we have now achieved a more homo- 
geneous coverage, especially filling the gaps in the radial range 
2'-5', in the area near the minor axis of NGC 4636. 

Figure[4] (middle panel) shows our estimate of the radial 
completeness of the kinematic GC sample: For radial bins of 1' 
width, we compute the ratio of GCs with velocity measurements 
to the number of candidate GCs for the published D+05 cata- 
logue (for consistency, GCs not listed in D+05 are not consid- 
ered here). Since bright GC candidates were preferred over faint 
ones for the spectroscopic observations, the completeness level 
changes significantly depending on the faint-end limiting mag- 
nitude, and the corresponding curves are shown with different 
symbols in the middle panel of Fig. [4] In the innermost bins, the 
completeness is very low, since the mask positions were chosen 
to avoid these parts where the light of NGC 4636 would domi- 
nate the GC spectra. As can already be seen from Fig. [1] (right 
panel), the spatial coverage in the area between about 2' and 7' 
is very good, and indeed the completeness peaks around a radial 
distance of 5.'5. Beyond about 8', the number of GCs with veloc- 
ity measurements becomes very small, hence the rapidly declin- 
ing completeness. Apart from the sparse spatial coverage, this 
is also due to the fact that the total number of GCs expected in 

2 The data set includes foreground stars which do not fulfil the colour 
criteria used for the GC candidate selection because all spectra, includ- 
ing those of the bright stars used for mask alignment, were extracted. 
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Fig. 3. Velocity distributions: Separating GCs from foreground stars. In both panels, the dashed line at 350 km s~' shows the velocity 
cut used in Paper I, and the grey region has a width of 100 kms 1 . In both panels, objects in this region are shown as white dots and 
numbered in decreasing order of the radial distance. The solid line indicates the systemic velocity of NGC4636 (v = 906 kms 1 ). 
Left: Colour vs. heliocentric velocity. The dashed box shows the range of parameters for the GCs analysed in Paper I (cf. Table 3 
therein). Large symbols show objects with velocity uncertainties Av < 65kms~'. The histograms (>'-axis labels in parenthesis) 
have a bin width of 50 kms -1 , and the thick line shows the objects with Av < 65 kms 1 . Right: Heliocentric velocity vs. projected 
galactocentric distance. Again, large symbols show objects with Av < 65 kms -1 . 
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Fig. 4. Spatial distribution of the velocity-confirmed GCs. Left panel: Radial distribution of the spectroscopic GC sample. The 
unfilled and grey histograms show the full sample and the data from Paper I, respectively. Middle panel: Radial completeness for 
different faint-end magnitude limits: For the GCs from the D+05 catalogue, we plot the ratio of GCs with velocity measurements 
with respect to the total number of GC candidates with colours in the range (0.73 < C—R < 2.15) found for the kinematic sample 
(cf. Fig. [5] left panel). Right panel: Azimuthal distribution of the GCs. The position angle (PA) is measured North over East and 
the dashed lines indicate the location of the photometric major axis of NGC4636. The histogram styles are the same as in the left 
panel. 



the outer regions is small, owing to the steeply declining number 
density profiles. This is also illustrated by Fig. |3](middle panel), 
where, for radii beyond 10' we find a number of foreground stars 
but almost no GCs. 



4.2. Luminosity distribution 

The left panel of Fig. [5] compares the luminosity distribution 
of the spectroscopic sample to the luminosity function (LF) of 
the GCs in the photometric catalogue. For a direct comparison 
between the two data sets, we only consider objects in the ra- 
dial range 2' < R < 8', i.e. where the spatial completeness 
of the spectroscopic sample is largest. For GCs brighter than 
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Fig. 5. Photometric properties of the GC sample. Left: Luminosity distribution of the GCs in the radial range 2' < R < 8' (indicated 
by the dashed lines in the middle panel of Fig. @). The grey histogram shows the distribution of the GC candidates from the D+05 
catalogue, the dashed bars are those D+05 GCs with velocity measurements. The unfilled histogram (thick black line) shows all 
GCs with 2' < R < 8' and Washington photometry and velocity measurements. The dashed line at = 23.33 mag indicates the 
turn-over magnitude (TOM) of the GC luminosity function (D+05). Middle: Colour-magnitude diagram. Small dots are the point 
sources from the D+05 catalogue (only objects within the radial range covered by our spectroscopic GC sample are shown), and 
open symbols are GCs with velocity measurements: Squares indicate those GCs with velocity measurements which are not in the 
D+05 list, circles are those identified in the D+05 catalogue. The solid line at C-R = 1.55 shows the colour adopted to divide red 
from blue GCs. The solid rectangle shows the area occupied by the velocity confirmed GCs in this study, and the vertical dashed 
lines at 0.9 and 2. 1 indicate the colour range adopted by D+05 to identify candidate GCs. Right: Colour distribution. The unfilled 
dashed histogram plots the colour distribution of the D+05 candidate GCs with magnitudes and radial distances in the range of our 
velocity-confirmed GCs. No statistical background has been subtracted; for graphical convenience, the number counts have been 
scaled by a factor of 0.5. The thick solid line shows the distribution for all velocity-confirmed GCs with Washington photometry, 
and the dashed histogram shows the colours of the velocity-confirmed GCs with R > 2.'5 (cf. Sect. l4.3l for details). 



niR ^21.4, the luminosity distributions are almost indistinguish- 
able indicating a high level of completeness, while, for fainter 
GCs, the sampling of the luminosity function becomes increas- 
ingly sparse. Note, that even the faintest GCs in the spectro- 
scopic sample are still ~ 0.6 mag brighter than the turn-over 
magnitude (wir tom = 23.33, D+05). 

4.3. Colour distribution 

The middle panel of Fig. [5] shows the colour-magnitude dia- 
gram of the NGC4636 GCS: Open symbols are the GCs from 
our spectroscopic study. Circles show GCs matched to the D+05 
final catalogue, and rectangles are the mostly bright GCs which 
are not in the D+05 catalogue (cf. Sect. 13. 11 1. Small dots show the 
objects from the D+05 catalogue (from the radial range covered 
by the spectroscopically confirmed GCs). The solid rectangle in- 
dicates the range of parameters of the GCs with velocities, and 
the dashed lines at C-R = 0.9 and 2.1 show the colour interval 
D+05 used to identify the photometric GC candidates. All but 
seven of our velocity-confirmed GCs have colours in the range 
used by D+05, confirming their choice of parameters. 

The right panel of Fig. [5] compares the colour histogram 
for the spectroscopic sample (thick solid line) to the distribu- 
tion of GC colours from the D+05 photometry (dashed unfilled 
histogram). Also for the smaller spectroscopic sample, the bi- 
modality is readily visible. The location of the peaks agrees well 
with the photometric GC candidates. For our sample, the heights 
of the blue and red peak are very similar, while, for the photo- 
metric sample, the blue peak is more prominent. This is due to 
the fact that the dashed histogram shows all D+05 photometric 
GC candidates for the full radial interval covered by our study: 
There is a stronger contribution from larger radii - where bl ue 



GCs dominate (see Sect. 110.41 ) - than for our kinematic sam- 
ple which becomes increasingly incomplete as one moves away 
from NGC 4636 (cf. Fig. S middle panel). 

4.4. GC colour distribution as a function of galactocentric 
radius 

Due to the different spatial distributions of metal-poor and 
metal-rich GCs (and the resulting differences with regard to the 
kinematics), a meaningful dynamical analysis requires a robust 
partition of the spectroscopic GC sample into blue and red GCs. 

D+05 used the minimum (C-R = 1.55) of the colour distri- 
bution of the GCs in the radial range 3 .'6 £ R £ S.'l to separate 
the two populations. 

We use the model-based mi xture modelling softwa r e pro- 
vided in the MCLUST package dFralev & Raftery||2002l 120061) 
to study the colour distribution of the NGC 4636 spectroscopic 
GC sample: MCLUST fits the sum of Gaussians (via maximum- 
likelihood) to the C-R data. The number of Gaussian components 
fit can either be specified by the user, or MCL UST co mputes the 
Bayesian Information Criterion (BIC, Schwarz [l978l) to find the 
optimal number of components. 

In Fig. [6] we plot the colours of the velocity-confirmed GCs 
versus their galactocentric distance (small dots). For a sliding 
window containing 105 GCs, we use MCLUST (in BIC mode) to 
determine the number of Gaussian components and plot the po- 
sitions of the peaks as function of radiufl (thin solid line). The 
large data points show the values obtained for the independent 

3 MCLUST is implemented in the R language and environment for sta- 
tistical computing and graphics (http: //www. r-pro ject.org). 

4 The derived quantities are plotted against the lower boundary of the 
radius interval covered by a bin. 
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Table 3. Colour distribution of the spectroscopic GC sample from normal mixture modelling 



ID 




N 


Model 


m R 


R 


N h 


Hb 


<Tb 


N T 




a, 


C-Rdivi 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


all 


458 


E 


19.11-2273 


0.51-15.43 


226 


1.25 


0.16 


232 


1.69 


0.16 


1.46 


2 


R > 2/5 


372 


E 


19.14-22.73 


2.53-15.43 


183 


1.26 


0.15 


189 


1.70 


0.15 


1.47 


3 


R < 2/5 


86 


X 


19.11-22.45 


0.51-2.50 


86 


1.42 


0.26 












R < 2/5 


86 


yt 


19.11-22.45 


0.51-2.50 


38 


1.22 


0.21 


48 


1.57 


0.19 


1.36 


4 


2.5 < R < 3/75 


86 


E 


19.64-22.62 


2.53-3.72 


49 


1.26 


0.16 


37 


1.72 


0.16 


1.50 



1 model type and number of components fixed. 

Note 3. The parameters were obtained using the MCLUST normal mixture modelling algorithm. The first Column describes and labels the subsets 
for which the analysis was performed. The total number of GCs is given in Col. 2, the type of model fit to the data is given in Col. 3, where 
'E' refers to equal-variance (homoscedastic) 2-component models, and 'X' to a single-component model. 'V is a heteroscedastic 2-component 
model. The range of i?-band magnitudes and radial distances are listed in Cols. 4 and 5, respectively. Columns 6 through 8 show the number of 
blue GCs, the position and width of the blue peak. Columns 9 though 1 1 are the same for the red GCs. Col. 12 is the colour dividing blue from red 
GCs. 
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Fig. 6. Colour distribution of the GCs as function of radius. The 
moving window contains 105 GCs, and the data points show the 
independent bins. The number of components to fit was deter- 
mined by the mclust software. The 'error-bars' give the widths 
of the respective distributions. The thick solid line shows the 
colour of the combined light of GCs within the sliding window. 
The vertical dashed line at 2/5 (^13 kpc) indicates the galacto- 
centric distance where the colour distribution becomes bimodal. 



bins, and the vertical bars indicate the width(s) of the distribu- 
tion^). The mean colour of the GCs is shown as thick solid line. 
Only for R ^ 2/5 (dashed vertical line), the colour distribution 
becomes bimodal. 

In Table[3] we list the parameters derived for various sub- 
samples of our data. Unless otherwise indicated, the number of 
components fit to the data was determined using the BIC. Apart 
from the GCs with R < 2/5 all distributions were found to be 
bimodal. 



Fig. 7. Colour distribution of the GCs within 2/5 (subset 3 in 
Table[3]). The thin solid curve is the Gaussian kernel density esti- 
mate for a bandwidth of 0.075 mag (same as the bin width of the 
histogram). The thick solid curve is the one-component model 
('X'), and the short-dashed curve is the unequal-variance 2— 
component model ('V'). The two components are shown with 
long-dashed curves, the corresponding peaks are indicated by 
the long-dashed vertical lines. For comparison, the thick grey 
line shows the kernel density estimate for the colour distribution 
of the GCs with R > 2/5 (number2 in Tabled, and the arrows 
mark the positions of the peaks. For reference, the colour used 
to separate blue from red GCs (C-R = 1.55, D+05) is shown as 
solid vertical line. 



4.5. Photometry of the GCs within 2! 5: 

As can be seen from Fig. [7] the colour distribution of the GCs 
with galactocentric distances smaller than 2/5 (sample 3 from 
Table[3) does not look bimodal, and MCLUST finds the distri- 
bution is best described by a single Gaussian (shown as thick 
solid line). Compared to the GCs with R > 2/5 (shown as thick 
grey curve), the distribution appears to be shifted towards the 
blue. When fitting a 2-component heteroscedastic model ('V in 
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Table[3) to the GCs within 2.'5, the fitted distributions (shown as 
thin dashed curves) are broader than for any of the other samples 
considered, and the peaks are offset to the blue. 

To test whether the MCLUST results for the GCs within 2. '5 
are due to the small sample size of 86 GCs, we apply the algo- 
rithm to the 86 GCs in the radial range 2 .'5 < R < 3 .'75 (sample 4 
in Table[3]l. For this sample of equal size, however, MCLUST de- 
termines that the most likely distribution is indeed bimodal with 
peak positions consistent with those found when considering all 
GCs outside 2.'5 (sample 2). 

We suspect that the photometry of the point sources within 
the central 2.5 arcmin is worse than that of sources outside that 
region because of issues with the subtraction of the galaxy light: 
The shape of the colour distributions is broadened by larger pho- 
tometric errors, and a small offset in the continuum subtraction 
results in a shift towards bluer colours. In the central 2.5 arcmin, 
it is therefore not possible to separate the blue and red GC pop- 
ulations. Consequently, the dynamical analysis of the subpop- 
ulations has to be restricted to galactocentric distances beyond 
2.'5(a 12.7kpc). 

Regarding the spatial distribution and the photometric prop- 
erties discussed above, we conclude that our spectroscopic sam- 
ple is a good representation of the GCs surrounding NGC 4636 
for galactocentric distances between 2.'5 and 8' (12.7 ^ R $ 
40.7 kpc). 

5. Definition of the subsamples 

As was demonstrated in Paper I, the presence of interlopers c an 
severely affect the derived line-of-sight velocity dispersion pro- 
file and, by consequence the inferred mass profile. In the follow- 
ing paragraphs, we describe our outlier rejection technique. The 
final subsamples to be used in the dynamical analysis are defined 
in Sect. |521 

5. 1 . Interloper rejection 

Objects with velocities that stand out in the velocity vs. galac- 
tocentric distance plot (Fig. [3] right panel) are potential outliers. 
Such 'deviant' velocities might be due to measurement errors, 
a statistical sam pling effect, or the presence of an intra-cl uster 
GC population iBergond et~aI1l2007t ISchuberfh et al.ll2008l) . In 
the low-velocity domain, possible confusion with Galactic fore- 
ground stars is the main source of uncertainty. 

5.1 .1 . Contamination by Galactic foreground stars 

In Sect. 13. 31 the division between foreground stars and bona-fide 
NGC 4636 GCs was made at v = 350kms~' (same as in Paper I). 
There are, however, seven objects with velocities between 300 
and 400 km s -1 (shown as labelled white dots in both panels of 
Fig. [3]). In order of decreasing distance from NGC 4636 these 
are: 

1. 7.1:15 (v = 313 ± 20^8"', m R = 19.83, C-R = 0.95) 
is almost certainly a foreground star: At a distance of 14.'8 
we hardly find any NGC 4636 GCs at all (cf. Fig. |3] middle 
panel). Moreover, this object would be very blue for a GC of 
this magnitude (cf. Fig. [5] middle panel). 

2. f09-57 (v = 345 ± 84kms"', C-R = 2.33, m R = 22.08) is 
probably a foreground star: Its very red colour lies outside 
the range of colours found for the GCs studied in Paper! 

3. Object f09-43 (v = 378 ± 45 km s _ \ m R = 20.82) has an ex- 
tremely blue colour (C-R = 0.78), and its velocity is offset 



from the NGC 4636 GC velocities found at this radial dis- 
tance (R = 10.'2), hence it is most likely a foreground star. 

4. f08-13 (v = 347 ± 66kms"', C-R = 0.95, R = 7.'9) is very 
blue and quite faint (m R = 22.5). 

5. 2.2:76 (v = 306 ± 60kms"', C-R = 1.40, m R = 21.9, 
R = 5.'3) is an ambiguous case. 

6. 1.2:15 (v = 392 ± 38kms _1 , C-R = 1.70, m R = 21.6, 
R = 4.'9) is an ambiguous case. 

7. f03-09 (v = 350 ± 69kms"\ C - R = 1.67, m R = 22.2, 
R = 4.'0) is also an ambiguous case. 

Thus, three of the seven objects in the velocity range 300- 
400 km s are almost certainly foreground stars because of their 
unlikely combination of extreme velocities and extreme colour. 
We remove these objects from the list of GCs. For the remaining 
four objects (No. 4-7) in this velocity interval, the situation is 
not as clear, so they remain in the samples to which we apply the 
outlier rejection scheme described below. 

5.1.2. Outlier rejection algorithm 

In this section, we apply the same outlier rejec tion method as 
descri bed in our study of the NGC 1399 GCS (ISchuberth et al.l 
2009) to our data. The method which is based on the tracer mass 
estimator bv lEvans et al.1 (120031) works as follows: For each sub- 
sample under consideration, we start by calculating the quantity 

1 N 

m N = -YjVt 2 ■ Ri , (1) 

/=i 

where v,- are the relative velocities of the GCs, and Rj their pro- 
jected galactocentric distances, and N is the number of GCs. 
Now, we remove the GC with the largest contribution to rtiN, 
i.e. max(v 2 • R) and calculate the quantity in Eq. [T]for the re- 
maining N - I GCs, and so on. In the upper panels of Fig. [8] 
we plot the difference between mj and m l+ \ against j, the index 
which labels the GCs in order of decreasing v 2 ■ R. 

For the blue GCs (shown in the right panel of Fig. [8} we con- 
sider only objects with velocity uncertainties Av < 65kms~'. 
The function mj - m.j+\ levels out after the removal of six GCs 
which are shown as dots in the lower panel of that figure. The 
thin solid curves enveloping the remaining GCs are of the form 

W(*)=^p, (2) 

where C max is the product v 2 • R for the first GC that is not re- 
jected. 

For the red GCs, the convergence of (mj - mj+\) is not as 
clear. We reject the two objects with low velocities which are 
most likely foreground stars (objects 6 and 7 from the list in 
Sect. 15.1.11 1. and find that the distribution of the remaining red 
GCs is symmetric with respect to the systemic velocity. 

5.2. The subsamples 

For the analysis of the kinematic properties of the NGC 4636 
GCS, we use the following four subsamples: 

- Blue All 209 blue GCs with R > 2.'5 

- BlueFinal 156 Blue GCs with R > 2.'5, Av < 65, six outliers 
removed 

- Red All 162 red GCs with R > 2. '5 
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Fig. 8. Outlier removal using the max(R ■ v 2 ) algorithm. Left: All red GCs. In the upper panel, we plot the quantity (ntj - nij+{) 
against the index j. The number of removed GCs ( j = 2) is indicated by the vertical dashed line and the arrow. In the lower sub- 
panel, we plot the heliocentric velocities vs. projected radius. Crosses mark GCs with velocity uncertainties Av > 65 kms -1 , and 
circles indicate GCs with Av < 65 km s~ . The rejected GCs are shown as dots. Right: The same algorithm applied to the Blue GCs 
with Av < 65km s . Six GCs are rejected. In both lower sub-panels, the systemic velocity of NGC4636 is indicated by a solid 
horizontal line. The vertical solid line at 2.'5 indicates the radius inside which the photometry does not permit a separation of blue 
and red GCs. The dotted vertical lines show the bins used for the dispersion profiles shown in Fig. Q2] 



- RedFinal 160 Red GCs with R > 2.'5, two outliers removed^ 6.1. The Anderson-Darling test for normality 

For reference, we compare these to the full sample (All) of when adopting p < 0.05 as criterion for rejecting the Null hy- 

459 GCs and the sample labelled AllFinal which combines pothesis of normality, we find that all subsamples are consistent 

BlueFinal and RedFinal and hence comprises 316 GCs. Details with bein 8 drawn from a normal distributionjThe^-values re- 

on the line-of-sight velocity distributions (LOSVDs) of these turned b Y the Anderson-Darling test (|Stephens||1974|) he in the 



samples are given in the following Section. 



range 0.15 < p < 0.72 (cf. Tableg]Col. 9). 



6. The line-of sight velocity distribution 

Figure [9] shows the line-of-sight velocity distributions 
(LOSVDs) for the samples defined above in Sect. 15.21 In 
each sub-panel, the unfilled histogram shows the respective 
initial sample, the filled histogram bars show the corresponding 
final samples. Below we comment on the statistical properties 
of these distributions which are compiled in Table|4] 



5 After removing the two outliers from the red GC sample, we find 
that further restricting the sample objects with velocity uncertainties 
below 60kms~' does not significantly change the velocity dispersion 
profile (cf. Table [6). Therefore, our final red sample still contains 13 
GCs with Av > 65km s~ 1 . 



6.2. The moments of the LOSVD 

For all samples, the median value of the radial velocities agrees 
well with the systemic velocity of NGC4636 (906 + 7kms~', 
Paper I) 

6.2.1. Velocity dispersion 

The velocity dispersion values quoted in Table d (Col. 6) 
were c alculated using the expressions given bv lPrvor & Meylan 
( 1993), in which the uncertainties of the individual velocity mea- 
surements are used as weights. 

The velocity dispersion of the final blue sample is 20 kms 1 
larger than that of the final red sample. Although this difference 
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Table 4. Statistical properties of the NGC 4636 globular cluster sample 



LLJ 




Ngc 


V 


V 


cr + Act 




K 


P \ i ^ iJ ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


All 




459 


916 


909 


197±7 


0.17±0.12 


0.05+0.23 


0.15 


AllFinal 


R > 2.' 5 & Av < 65 & OutRej 


387 


920 


912 


181+7 


0.11+0.09 


-0.45+0.13 


0.15 


Blue 


C-R < 1.55 &«> 2.'5 


209 


928 


912 


198±11 


0.30±0.15 


0.03+0.28 


0.31 


BlueFinal 


C-R <\.55&R> 2.'5 & Av < 65 & OutRej 


156 


925 


912 


173±10 


-0.10+0.12 


-0.69+0.18 


0.17 


Red 


C-R> 1.55&«>2.'5 


162 


893 


900 


162+10 


-0.03+0.20 


0.12+0.40 


0.64 


RedFinal 


C-R > 1.55 &R > 2.'5 & OutRej 


160 


900 


901 


153+10 


-0.03+0.14 


-0.48+0.20 


0.72 



Note 4. Cols. (1) and (2) label the samples, and the number of GCs is given in (3). Cols. (4) and (5) list th e mean and median 
radial velocity (in kms"'). The velocity dispersion in (kms -1 ) as computed with the Prvor & Mevlan (1993) formula is given in 
Col. (6). Cols. (7) and (8) give the skewness and the reduced kurtosis, the uncertainties were estimated using a bootstrap with 999 
resamplings. The p-value returned by the Anderson-Darling test for normality is given in Col. (9). 
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Fig. 9. NGC 4636 line-of-sight velocity distributions. From top 
to bottom, the panels show the velocity distribution for the entire 
sample, the blue, and the red subsample, as defined in Sect. 15.21 
In all panels, the vertical dashed line shows the systemic veloc- 
ity of NGC 4636. The unfilled histograms show the respective 
velocity distributions prior to any quality selection and outlier 
rejection. The filled histograms show the final samples. In all 
panels, the dashed curve shows the Gaussian corresponding to 
the dispersion value for the respective final sample as given in 
Table|4] The thin solid line is the corresponding Gaussian kernel 
density estimate (for a bandwidth of 50 km s , same as the bin 
width). 



butions are symmetric with respect to the systemic velocity of 
NGC 4636. The only (significantly) non-zero skewness is found 
for the blue GCs prior to weeding out outliers and GCs with 
large measurement uncertainties. 

6.2.3. Kurtosis 

Column 8 of Table|4] lists the reduced kurtosis (i.e. a Gaussian 
distribution has k — 0) of the respective subsamples. The fourth 
moment of the LOSVD reacts quite severely to the treatment 
of extreme velocities: The removal of only two clusters from 
the Red sample changes the kurtosis from k = 0.12 ± 0.40 to 
k = -0.48 + 0.20 (RedFinal). The FinalBlue sample also has 
a negative kurtosis, meaning that both distributions are more 
'flat-topped' than a Gaussian. Indeed, the kernel density esti- 
mates (thin solid lines in Fig. [9]) have somewhat broader wings 
and a flatter peak than the corresponding Gaussians (thick dotted 
curves). These differences, however, are quite subtle and we sus- 
pect that the samples considered in this work are probably still 
too small to robustly determine the 4 th moment of the velocity 
distributions. The slightly negative kurtosis values are, however, 
consistent with isotropy (a projected kurtosis of zero is expected 
only for the isothermal sphere). 

We conclude that the final GC samples that will be used for 
the dynamical analysis are symmetric with respect to the sys- 
temic velocity of NGC 4636 and do not show any significant de- 
viations from normality. 

7. Rotation 

Due to the very inhomogeneous angular coverage of the data, 
the results from the search for rotation of the NGC 4636 GCS 
presented in Paper I were quite uncertain. Below, we use our en- 
larged data set for a re-analysis and compare the findings to the 
ones in Paper! To detect signs of rotation, we fit the following 
relation to the data: 



cannot be considered significant (since the values derived for 
blue and red GCs marginally agree within their uncertainties), it 
will be shown below (Sect. [8}, that the radial velocity dispersion 
profiles of the two subpopulations, however, are very different 
(cf.Figs.[T2]and[]3D. 

6.2.2. Skewness 

The skewness values (Table|4] Col. 7) for the final blue and red 
samples are consistent with being zero, i.e. the velocity distri- 



v,(®) = v sys +Asin(@-@ ) , (3) 

where \\ is the measured radial velocity at the azimuth angle 0, 
A is the amplitude (in units of kms ), and Qp i s the position 
angle of the axis of rotation (see Cote et al.ll200ll for a detailed 
discussion of the method). 

In Table|5J we present the results of the rotation analysis for 
the different subsamples of our data: Columns 2 and 3 give the 
values of ©o and the amplitude A found for the samples pre- 
sented in Table|U Columns 4 to 6 give the corresponding results 
obtained for the GCs in the radial range 2.'5 < R < 8.'0, where 
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Table 5. Rotation of the NGC 4636 globular cluster system 











2.'5 < R < 8' 




Sample ID 


®o 


A 


Ngc 


©0 


A 




[deg] 


[kms" 1 ] 




[deg] 


[kms" 1 ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


All 


160+21 


34±14 


325 


152+16 


49+16 


All. Final 


150±15 


46±15 


269 


141+12 


59+15 


Blue 


156±18 


58+21 


185 


144+17 


66+23 


Blue. Final 


144+13 


77±21 


141 


140+12 


88+22 


Red 


172±27 


37+19 


138 


175+35 


32+20 


Red. Final 


171+26 


36+18 


136 


174+34 


30+19 



Note 5. The position angle of NGC 4636 is 150°. The first column gives 
the sample identifier (Sect. l5T2t . Columns 2 and 3 give the axis of ro- 
tation and the amplitude, respectively. Columns 4 through 6 give the 
number of GCs, rotation angle and amplitude for the samples restricted 
to radii R < 8' (i.e. the range in which the completeness of the spectro- 
scopic GC sample is large). 



Table 7. 

bins 



GC Velocity dispersion profiles for constant number 



Blue GCs (BlueFinat) 



No. 


R 


cr±Acr 


n 


R 


cr±Ao- 


n 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




1.'20 


208+30 


26 










2.'03 


242+35 


26 


1/85 


235+31 


30 


1 


3.'05 


202+24 


39 


3.' 13 


135+16 


40 


2 


4.' 55 


179+21 


39 


4/44 


152+19 


40 


3 


5/75 


162+20 


39 


5/81 


172+21 


40 


4 


7/92 


143+19 


39 


8/48 


151+19 


40 



Red GCs (RedFinal) 



Note 7. The first column numbers the bins (for GCs with R > 2/5). 
Column (2) gives the mean galactocentric distance of the blue GCs, 
Col. 3 gives the velocity dispersion in units of kms~', and Col. 4 is 
the number of blue GCs in that bin. Columns 5 through 8 give the cor- 
responding values for the red GCs. For completeness, the first two rows 
give the values for GCs inside 2/5 where the colour distribution is not 
bimodal. 



the spatial completeness of our sample is highest (cf. Sect l4. 11 1. 
It appears that the rotation signal (within the uncertainties) is 
robust with respect to the radial range considered and the appli- 
cation of the outlier rejection algorithm. 

To search for variations of the rotation signal with galacto- 
centric radius, we plot, in the left and middle panel of Fig. [TO] 
the rotation parameters ©o and A determined for moving bins 
of 50 GCs. The following paragraphs summarise the results for 
blue and red GCs. 

7.1. Rotation of the blue GCs 

Table|5]shows that the rotation signature is strongest for the blue 
GCs (final sample). 

In the right panel of Fig. [10] we plot the radial velocities 
against the position angle for the final blue sample restricted to 
galactocentric distances below 8' (which, in Table[5]is the sam- 
ple with the strongest rotation signature). The data for the 141 
GCs are shown as circles, and the squares mark the mean veloc- 
ity calculated for 60° wide bins. The thick solid line shows Eq.[3] 
(A = 88kms~' and O = 140°). 

7.2. Rotation of the red GCs 

For the red GC samples, the overall rotation signature as quoted 
in Table[5]is weaker than that of the blue GCs. Within R ^ 6/0, 
the amplitude for the final sample plotted in the middle panel 
of Fig. [10]is consistent with being zero (104 GCs with 2/05 55 
R & 6/0 : A = 17 ± 22kms" 1 ,0 o = 186 + 73°). Only the 
last few bins suggest an increase with the amplitude reaching 
values of ~ 80kms"' (32 GCs with 6/0 < R < 8/0 : A = 
82 ± 37km s" 1 , 0o = 160 ±21°). 

Between 3'and 5', the axis of rotation changes from being 
aligned with the minor axis to the photometric major axis of 
NGC 4636. Given the low values of the amplitude at these radii, 
however, this change of the axis of rotation remains uncertain. 
For radii beyond 5', the axis of rotation remains constant and 
coincides with the major axis. 

7.3. Comparison to Paper I 

The findings in this section deviate significantly from the values 
presented in Paper I, where no rotation was detected for the blue 



GCs, while the strongest signal was found for the red clusters 
within 4' (©o = 72 + 25°, A = -144 + 4410ns" 1 ). 

The discrepant findings are likely due to the azimuthal in- 
completeness of the data in Paper I. 

8. Globular cluster velocity dispersion profiles 

FigureQ~2]shows the velocity dispersion profiles obtained for the 
blue (upper panel) and red (lower panel) subsamples defined in 
Sect. 15.21 The GC velocities were grouped into six radial bins 
(shown as dotted lines in the lower panels of Fig.[8j. The first bin 
comprises the GCs within 2/5, and the following four bins (start- 
ing at 2/5, 4/0, 5/5 and 7/0) have a width of 1 .'5. The outermost 
bin (8/5 < R < 15/5) collects the GCs in the more sparsely pop- 
ulated (and sampled) outer GCS. The data are given in Table|6] 
Circles are the values obtained for the Blue and Red subsamples 
and dots represent the dispersion profiles determined for the fi- 
nal samples {BlueFinal and RedFinal). The values from PapeiT 
(Table 4 therein) are shown as diamonds. 

8.1. Blue GCs 

Compared to the initial sample, the removal of GCs with veloc- 
ity uncertainties Av > 65 km s _1 and the six probable interlopers 
identified in Sect. 15. 1 .21 significantly reduces the velocity disper- 
sion, in particular in the 4th and 5th bin. The values for the final 
blue sample (dots) agree very well with the data from Paper I 
(diamonds). The dispersion profile of the blue sample declines 
with galactocentric radius. 

8.2. Red GCs 

The velocity dispersion profile of the red GCs (cf. Fig. [12] 
lower panel) shows a more complex behaviour. The sudden drop 
by almost 100 kms~' from the first to the second bin was al- 
ready present in the data for Paper I. While the poorer data in 
Paper I still left the possibility of a sampling effect, our new 
enlarged data shows the same feature which we attribute to the 
fact that blue and red GCs cannot be separated within the cen- 
tral 2/5: Inside this radius, the presence of a substancial num- 
ber of 'contaminating 'metal-poor GCs (which would have a 
broader velocity distribution) within our 'red' sample dominates 
the measurement and leads to the very high dispersion value. 



12 



Y. Schuberth et al.: Dynamics of the NGC4636 globular cluster system II 



360 
270 

a) 

S180 
© 

90 


200 

T |100 

< o 

-100 





iii J 


L 


Blue GCs I 






L 




= 


-i%f* 1 


1 — 1 — h 






H — 1 — — 1 — 1 — | — 1 — — 1 — f 








mm m i 











360 
270 



2180 
© 



90 -V- 



Red GQs 




200 

g 100 

< 

-100 



1500 : 



1200 



S900 



>600 



300 



10 



10 



R [arcminj 



R [arcminj 




180 



Fig. 10. Rotation of the NGC 4636 globular cluster system. Left: Circles show the rotation of the final blue sample {BlueFinal) 
computed for moving bins of 50 GCs with a step site of 5 GCs. The grey area indicates the radial coverage and the uncertainties. 
Dots indicate independent bins. The upper sub-panel shows the rotation angle 0o vs. projected galactocentric radius. The NGC 4636 
photometric major and minor axis are shown as long-dashed and dotted horizontal lines, respectively. The lower sub-panel shows 
the amplitude A as function of radius. Middle: The same as the left panel but for the final red sample (RedFinal). In both plots, 
the dashed vertical lines indicate the radial range 2. '5 < R < 8' for which we have the best spatial coverage of the GCS (see 
also Sect: l4. lb Right: Velocity versus position angle for the 141 blue GCs (circles) from the sample BlueFinal with distances 
between 2.' 5 and 8'. The thick solid curve shows the best fit of Eq.[3]to the data, the thin dashed curves are offset by the sample 
dispersion (cr = 180 kms 1 ). The vertical dashed and solid lines show the best-fit ©o and the photometric major axis of NGC 4636, 
respectively. The grey areas indicate the uncertainty AQq. The squares show the mean velocity for bins of 60°. The unfilled histogram 
at the bottom shows the angular distribution of the GCs (the bins have a width of 20° and the counts are multiplied by a factor of 20 
for graphic convenience). 



Table 6. Velocity dispersion profiles for fixed radial bins 



Blue 



BlueFinal 



Red 



RedFinal 



N° 


range 


ii 


<r 


n 


(T 


n 


(T 


n 


cr 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 


O.'O- 2.'5 


56 


249±25 


52 


226±23 


30 


235+32 


30 


235+32 


2 


2.'5- 4.'0 


55 


199+21 


44 


195+22 


45 


140+16 


45 


140+45 


3 


4.'0- 5.'5 


60 


200±20 


45 


186+21 


50 


183+21 


48 


154+48 


4 


5.'5- 7.'0 


51 


200±23 


38 


146+21 


32 


169+23 


32 


169+32 


5 


7.'0- 8.'5 


25 


206±32 


18 


154+28 


22 


165+27 


22 


165+22 


6 


8.'5- 15/5 


18 


166+33 


11 


131+34 


13 


122+31 


13 


122+13 



Note 6. Column 1 gives the bin number and Col. 2 is the radial range covered by a given bin. Col. 3 is the number of GCs in the Blue subsample, 
Col. 4 the lineof sight velocity dispersion in kms~'. Cols. 5 and 6 are the same for the BlueFinal sample. The corresponding values for the red 
GCs are given in Cols. 7 through 10. Columns 1 1 and 12 show the values obtained when further restricting the RedFinal sample to velocities with 
uncertainties below 3 < 65km s . 



To search for colour trends amongst the GCs inside 2. '5, we 
plot in FigfTTKupper panel), the the line-of-sight velocity disper- 
sion as a function of colour. Dots represent the 82 GCs within 
2.'5 with velocity uncertainties below 65 kms _I . The 316 GCs 
from the Final sample are shown as squares. Inside 2.'5, there 
is no discernible trend with colour, except for a dip near the 
colour used to divide blue from red GCs. For GCs outside 2. '5, 
we observe a constant dispersion for the red GCs, which then 
increases towards bluer colours. However, there is no sign of a 
' jump' near the dividin g colour as seen in the NGC 1399 GCS 
(Schubert h et al.l [20101) . For GCs with projected galactocentric 
distances between 2.'5 and 4.'0, we find a low dispersion of only 
140 ± 16 km s . This value then again rises to almost 170 km s 
between 5. '5 and 8.'0. 



9. NGC 4636 stellar kinematics 

Long slit sp ectra of NGC 4636 were obtained by iBender et al.l 
(119941) . and iKronawitter et al.l d2000t) p resented d e tailed mod- 
elling based on these data. Recently, iPu & Hani d201 ll) pre- 
sented radial velocities and velocity dispersions derived from 
deep long-slit spectra along the major and minor axis. Their data 
reach out to ~ 2.'3, corresponding to about 1 1.5 kpc. 

We will use the stellar kinematics to constrain the halo mod- 
els derive d for the GCs (see Sect. 111.1.3) . Below, in Fig. [13] we 
show the lBender et al.l (U994) data and compare them to our own 
measurements obtained in parallel with the GC observations. 

9. 1 . FORS 2 spectra of NGC 4636 

During the first MXU observations of the NGC 4636 GCs, we 
placed 29 slits (along the North-South direction) on NGC 4636 
itself. This mask (Mask 1_1 from Paper I) has a total exposure 
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Fig. 11. Upper panel: Line-of-sight velocity dispersion as a 
function of colour. Dots represent the values for the 82 GCs 
within 2.'5 and velocity uncertainties below 65 kms -1 , squares 
show the results for the 316 GCs of the AllFinal sample. The la- 
bels give the number of GCs per bin. The bin size is 0.15 mag, 
but the bluest and reddest bins have been resized so that no ob- 
jects are excluded. The horizontal bars mark the colour range 
covered by a given bin. Lower panel: Heliocentric velocity as 
a function of colour. The symbols are the same as in the upper 
panel. 



ctq = 243 + 3kms (shown as diamond in the right panel of 
Fig. [T3l . The low central dispersion quoted by BSG94 is likely 
due to the instrumental setup: These authors used a slit of 2'.' 1 
width, and at this spatial resolution the luminosity weighted dis- 
persion measured for the centre may be subst antially lower than 
values obtained using a smaller slit width: iProctor & Sansorrj 
used a slit of l'/25 which provides a similar spatial resolution 
as our l'.'O wide MXU slits, yielding similar dispersion values. 

For the dynamical modelling, dispersion values at large ra- 
dial distances are of particular interest. Unfortunately, due to the 
low S/N in the small slits we used, the quality of our data de- 
grades for radial distances beyond R ^ 0.'7: The uncertainty 
of the individual data points increases and so does the scatter. 
However, the velocity dispersion seems to decline as indicated 
by the solid curve in Fig. [13] (right panel) which shows a linear 
fit to <t 1os (7?) for R > 0.'5. T his trend is confir med by the recently 
published measurements bv lPu & Hani (1201 ll) (small dots in Fig. 
[T3l which clearly show a declining stellar velocity dispersion for 
radii beyond ~ 0.'6 (~ 3 kpc). These data have a higher S/N than 
our FORS measurements and extend out to 2.'3, i.e. almost into 
the regime where GC dynamics is available for both the metal- 
poor and the metal-rich subpopulation (shown as filled triangles 
in Fig. [T31 ). In this context, it is interesting to note that the very 
low velocity dispersion of 135 ± 16 km s , observed for the red 
GCs near ~3.'l is consistent with the outermost stellar veloc- 
ity dispersion value of 138 ± 26 kms -1 . This might suggest a 
connection between star s and metal-rich GCs similar to the one 
reported for NGC 1399 lSchuberth et al.ld2010b . 

For the dynamical modelling of the st ellar component of 
NGC 4636 we use the iBender et ail d 19941) data in the radial 
range 0.'05 to 0.'6 (0.25-3.1 kpc), indicated by the dashed lines 
in the right panel of Fig. [13] The central data points are not in- 
cluded in our modelling since the deprojection of the luminosity 
profile and, by consequence, the stellar mass profile are only re- 
liable for R > 100 pc * 0.'2 (see Sect. [1021 . 



time of 2 hours (7 200 s), and the final image is the co-addition 
of three consecutive exposures, thus ensuring a good cosmic-ray 
rejection. The data were reduced in the same manner as the slits 
targeting GCs. The sky was estimated from the combination of 
several sky-slits located ~ 3' from the centre of NGC 4636. 

Using the pPXF (penali sed PiXel Fitting) routine by 
Cappellari & Emsellem (2004), we determined the line-of-sight 
velocity dispersion for the one-dimensional spectra. The uncer- 
tainties were estimated via Monte Carlo simulations in which 
we added noise to our spectra and performed the analysis in the 
same way as with the origina l data. The templates used in the 
analysis were taken from the IVazdekisI (l999) library of syn- 
thetic spectra. Our data are shown in Fig. [13] where unfilled 
and filled squares refer to slits placed to the North and South 
of the centre of NGC 4636, respectively. The data are listed in 
TablelAJl 

9.2. The stellar velocity dispersion profile of NGC 4636 

In Fig. [13] we comp are our results (squares) with the values from 
Bend er et all (1 19941) (large dots). Within the radial range 0.'04 ^ 
R 5 0.'6, the agreement between both data sets is excellent. 

For the central velocity dispersion (measured from a 2" 
long slit within about 1" of the centre of NGC 4636), we find 
cr = 233 ± 6kms~'. This is substantially higher than the 
value published by BSG94 (cr = 211 ± 7kms _1 ). However, 
the high central velocity d ispersion we find is suppo rted by the 
measurement published bv lProctor & Sansom (2002) who found 



10. Jeans models for NGC 4636 

In the next paragraphs, we give the relevant analytical expres- 
sions and outline how we construct the spherical, non-rotating 

Jeans models for NGC 4636. 

In PapeiT, we chose an NFW-halo (iNavarro et al.lll997h to 
represent the dark matter in NGC 4636. In this work, we will 
consider both NFW halos and two mass distribut ion with a finite 
central density: The cored profile proposed by Burkert (1995) 
which has the same asymptotic behaviour as the NFW halo and 
the logarithmic potential which leads to (asymptotically) flat ro- 
tation curves. 

10.1. The Jeans equation and the line-of-sight velocity 
dispersion 

The spherical, non-rota ting Jeans equation (see 
e.g. lBinnev & Tremaindll987l) reads: 

d(n(r)<7») B(r) 2 G- M(r) 

j + 2 n(r) cr r (r) = -n(r) , (4) 

dr r r l 

erg 2 

with B = 1 - -V . 

°7 

Here, r is the radial distance from the centre and n is the spatial 
(i.e., three-dimensional) density of the GCs; o> and cr g are the 
radial and azimuthal velocity dispersions, respectively. /3 is the 
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anisotropy parameter, M(r) the enclosed mass (i.e. the sum of 
stellar and dark matter) and G is the constant of gravitation. 

For our analysis , we use the expressions given by, 
e.g. Mamon & Lokas (2005), see also van der Marel & Franx 
dl993|). Given a mass distribution M(r), a three-dimensional 
number density of a tracer population n{r), and a constant 
anisotropy parameter/?, the solution to the Jeans equation (Eq.|4]l 
reads: 

n(r)oj.(r) = G J n(i)M(j)-|-j ds. (5) 
This expression is then projected using the following integral: 



N{R) 



J R Vr 2 - R 2 Jr 



finer 1 , dr 



r^r 2 -R 2 



(6) 



where N(R) is the projected number density of the tracer pop- 
ulation, and cr los is the line-of-sight velocity dispersion, to be 
compared to our observed values. In the following, we discuss 
the quantities required to determine cri os (R). 



where T* « is the 7?-band mass-to-light ratio of the stellar pop- 
ulation of this galaxy (see below in Sect. 110.3b . 

Since the integral in Eq. [10] cannot be expressed in terms 
of simple standard functions, we use an approximation in our 
calculations: The inner part (r ^ 45 kpc) is represented by a 
sequence of polynomials, while the behaviour at larger radii is 
well represented by an arctan function. The stellar mass profile is 
plotted in Fig. [15] and the expressions and coefficients are given 
in AppendixlBl 

10.3. The stellar mass-to-light ratio 

In Paperl, we used an /?-band Y* = 6.8. This was de- 
rived from the dyn a mical estimate for the B-band given 
by iKronawitter et al.l (I2000'). Having adopted a distance of 
17.5 Mpc for our current analysis, the value from Paperl which 
was based on a distance of 15 Mpc is reduced to = 5.8. 
We will adopt this value for our dynamical modelling of the 
NGC 4636 GCs. 



10.2. Luminous matter 

To assess the stellar mass of NGC 4636 we need to deproject 
the galaxy's surface brightness profile. Moreover, to consistently 
model the line-of-sight velocity dispersion profile of the stars (cf. 
FigfTSl. we require analytical expressions for both, the projected 
and the three-dimensional stellar density. 

As in Paperl, we use the data published by D+05 (shown as 
dots in the upper panel of Fig. [T4b . for which the authors gave 
the following fit: 

//(/?) = 

-2.5 log (3.3 ■ 10- 7 (1 + + 5.5 ■ HT* (l + £f 13 ) . 

(7) 

Their fit is shown as dashed line in Fig. [14] There is, however, 
no analytical solution to the deprojection integral for this func- 
tion. We therefore fit the data using the sum of three Hubble- 
Reynolds profiles instead: 



( 3 



M#) = -2.51og 



V/=l 



1 



R 

Ro,i 



(8) 



where the parameters are given in Table[8] Our fit is shown as 
solid black line in Fig. [14] and the thin gray lines indicate the 
three components. The deprojection of Eq. [8]reads: 



j(r) 



pc J 



= y Wo 



1 + 



R', 



-(or, + 1/2) 



(9) 



Where S is the Beta function and N'qj = Cm r -Nqj, where Cm r 
2.1') 2 x 10 10 is the factor converting the surface brightness into 
units of L pc~ 2 for M q r = 4.28. The radii are in pc, i.e. for a 
distance of 17.5 Mpc R' QJ = 5.09 x 10 3 ■ R QJ . 
The lower panel of Fig. [14] compares the deprojection as given 
in Eq. [9] (solid black line) to the curve obtained by numerically 
deprojecting Eq.[7] Within the radius interval covered by the data 
points, both deprojections agree extremely well, and we proceed 
to use the analytical expression given in Eq. [9] to represent the 
density distribution of the stars in NGC 4636. 

The stellar mass profile is then obtained through integration: 



M(r) = Y*, R ■ 4tt j(s)s 2 ds 



r 

Jo 



(10) 



10.4. Globular cluster number density profiles 

Below we present the fits to the number density profiles of the 
GC subpopulations as listed in Table A3 of D+05 and the ana- 
lytical expressio ns for the deprojection s. As in our study of the 
NGC 1399 GCS dSchuberth et al .12 010) we parametrise the two- 
dimensional number density profiles in terms of a Reynolds- 
Hubble law: 



N(R) = N 



-5 



(11) 



where Ro is the core radius, and 2 • a is the slope of the power- 
law in the outer region. For the above expression, the Abel inver- 
sion has an analytical solution and the three-dimensional num- 
ber density profile reads: 



(r) = 



No 



1 



Ro 



1 + 1- 

Ro 



2\ 



-(*+!) 



(12) 



where S is the Beta function. For both subpopulations, the fits 
are performed for the radial range 2. '5 < R < 13. '5 where the 
lower boundary is the minimum radius where blue and red GCs 
can be separated. The upper boundary corresponds to the ra- 
dius where the GC counts reach the background level (D+05). 
The parameters obtained for the blue and red GCs are given in 
Table[8] Figure[l6] shows the data and the fitted profiles. Note 
that the profile of the red GCs is significantly steeper than that 
of the blue GCs. 

10.5. The dark matter halo 

All three dark matter halos considered in this study have two free 
parameters, allowing a direct comparison of the results. 

10.5.1. The NFW profile 

The mass profile of the NFW halo reads: 



M NFW (r) = 4ng s r s 



(13) 



where g s and r s are the characteristic density and scale radius, 
respectively. 
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Table 8. Fit parameters for the luminosity density profile of 
NGC 4636 and the GC number density profiles 





N 




a 


*(§.«) 


Lumprof 1 


2.07 x irr' 


3.47 x 1CT" 


1.1 


1.887 


Lumprof 2 


9.73 x lfr 8 


1.47 x 10-' 


1.1 


1.887 


Lumprof 3 


5.45 x lfr 9 


1.41 


1.3 


1.708 


Blue 
Red 


12.0±1.1 
12.4+1.1 


6.0±0.3 
6.8±0.3 


2.0±0.1 
3.0+0.1 


4 
3 
16 
15 



Note 8. The first three rows give the parameters for the three- 
component Hubble-Reynolds profile fit to the luminosity density profile 
(Eq.[8). The last two rows are the parameters found for the GCs (Eq.lllt 



To express the halo parameters in t erms of concentration and 
virial mass, we use the definitions from Bullock et al. (2001) and 
define the virial radius R v [ r such that the mean density within 
this radius is A v ; r = 337 times the mean (matter) density of the 
universe (i.e. 0.3 p c ), and the concentration parameter is defined 
as c v ir = R Ylt /r s . 

10.5.2. The Burkert halo 



The density profile for the cored halo which iBurkertl (ll995l) in- 
troduced (to represent the dark matter halo of dwarf galaxies) 
reads: 

Q(r) = r , (14) 

c + *)H) 

and the cumulative mass is given by the following expression: 
M(r) = 







♦lb 
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r 2 ' 


1 / 










1 + 




arctan 












r 2 


2 \ 





(15) 



10.5.3. The logarithmic potential 

The lo garithmic potential (see lBinnevll98UlBinnev & Tremairiel 
Il987l) . by construction, yields (asymptotically) flat rotation 
curves. In contrast to the NFW profile, it has a finite central den- 
sity. The spherical logarithmic halo has two free parameters, the 
asymptotic circular velocity Vo, and a core radius Tq. The mass 
profile reads: 

M Log (r) = - • -2- , (16) 

where G is the constant of gravitation. 



10.6. Modelling the velocity dispersion profiles 

To find the parameters that best describe the obse rved GC veloc- 
ity dis persion data, we proceed as described in Schu berth et al.l 
(2009): For a given tracer population and anisotropy B e 
{-0.5, 0, +0.5), we create a grid of models where the density (or 
vo in case of the logarithmic halo) acts as free parameter while 
the radii have discrete values, i.e. r<jark e {1,2,3, 100} kpc. 
For each point of this grid, the line-of-sight velocity dis- 
persion (Eq. |6]l i s com puted using the expressions given in 
M amon & Lokasl (l2005h . where the upper limit of the integral 
in Eq.|6]is set to 600 kpc. 



To find the joint solution for the different tracer populations 
(labelled a and b), we determine the combined parameters by 
minimising the sum;^ 2 = xl + x\- The confi dence level (CL) 
contours are calculated using the definition bv lAvnil (1 19761) . i.e. 
using the difference 5x 2 above the minimum^ 2 value. With two 
free parameters, e.g. (r,jark>Pdark) the 68, 90, and 99 per cent con- 
tours correspond to x 2 = 2.30, 4.61, and 9.21, respectively. The 
results for the three dark matter halos (NFW, Burkert and loga- 
rithmic potential) are presented in the following section. 



1 1 . The mass profile of NGC 4636 

We model the observed line-of-sight dispersions for the final 
samples of the red and blue GCs shown in Fig. [19] and listed in 
TableQfor three different parametrisations of the dark halo. 



11.1. Jeans models for an NFW halo 

11.1.1. NFW Halo: Results for the blue GCs 

In Fig. [17] we show the NFW models for the blue GCs after 
transforming the parameters to th e (M v ; r , c v ; r ) plane using the 
definitions in Bullo ck et alJ d200l[) . In all panels, the respective 
best fit value as given in Table[9]is shown as a cross. Circles are 
the corresponding values from Paper I. In all cases, these values 
lie within the 68% CL contour of the present study. 

From the very elongated shape of the confidence level con- 
tours it is apparent that while the GCs can be used to estimate 
the total mass of the halo, the concentration is only poorly con- 
strained. As will be shown in Sect. II 1.4~T1 this degeneracy can be 
partially overcome by considering models for the stellar velocity 
dispersion profile of NGC 4636. 

The best-fit dispersion profiles for the three values of the 
anisotropy parameters (p 1 e {-0.5, 0, +0.5}) are shown as thin 
black lines in the left panel of Fig. [19] and the corresponding 
parameters of the NFW halos are listed in Table|9] (Cols. 4— 
9). For all three values of B, a very good agreement between 
data and models can be achieved, and the differences between 
the x 2 values is marginal. The models diverge at small radii 
(R £ 2.' 5 ~ 13 kpc) where the velocity dispersion and the shape 
of the number density profile of the blue GCs cannot be well 
constrained. For comparison, we plot (as dot-dashed line) the 
velocity dispersion curve expected if there were no dark matter 
and the only mass were that of the stars. 

As expected, the best-fit halo derived assuming a tangential 
orbital anisotropy (B = -0.5, B.tan, long-dashed line) is less 
massive than the one obtained in the isotropic case (B.iso, solid 
line) and the model for a radial bias (B = +0.5, B.rad, short- 
dashed line) returns the most massive dark halo. This is also 
illustrated by the bottom left panel of Fig [19] where the corre- 
sponding mass profiles (thin black lines) are shown in terms of 
the circular velocity. 

1 1 .1 .2. NFW halo: Results for the red GCs: 

The resulting halo parameters for the red clusters are listed in 
Table[9](Cols. 4-9) and illustrated in the middle panel of Fig. [19] 
Here, the agreement between data and models is worse than in 
the case of the blue GCs. A considerable part of the uncertainty 
is caused by the curiously low value at ~16 kpc. In spite of this, 
the resulting circular velocities (shown as thick grey lines in the 
bottom left of Fig. [T9l of the different halo models are not dra- 
matically different from those of the blue GCs. For all three 
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Blue GCs,/? = -0.5 Blue GCs, /J = Blue GCs,/? = +0.5 
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M^tMJ fcUfMJ M^IMJ 



Fig. 17. Jeans models for the blue GCs (BlueFinal) for an NFW-type dark halo. From left to right, the panels show the results for 
B = -0.5,0 and +0.5. The parameters are shown in the (M v ; r , c v ; r )-plane. The thick solid, dashed and thin solid lines indicate the 
68,90, and 99 per cent confidence limits. The colour map is the same for all panels, and the cross indicates the location of the 
minimum^ 2 value. The square shows the minimum^ 2 for the joint models obtained from the blue GCs and an isotropic model for 
the NGC4636 stellar dispersion profile. The circle indicates the respective best-fit value from Paper! All model parameters are 
listed in Table[9] In all panels, th e long-dashed (dash-dotted) lines show the median (68 per cent values) for simulated NFW halos 
as found by Bulloc k et all (l200ll) . 



values of B, the circular velocity stays approximately constant 
within 40 kpc. 

11.1.3. Model for the stars 

W e use the stellar velo city dispersion measurements presented 
bv lBender et alJ (1 19941) to constrain the concentration parameter 
of the NFW halo. The right panel of Fig.|20]shows the (M vir -c vir ) 
plane for the isotropic case. High concentrations are excluded, 
since adding large amounts of dark matter in the central parts of 
NGC 4636 would severely overestimate the velocity dispersion 
profile of the stellar component. 

1 1.2. Jeans models for a Burkert halo 

Figure[18]shows the parameter space explored to find the best-fit 
isotropic Jeans model for the blue GCs for a Burkert-type dark 
halo. The best-fit Burkert models for the GCs are shown in the 
middle panels of Fig. [19] and the parameters are given in Cols. 
11-12 of Table [9] The circular velocities corresponding to the 
different mass distributions are compared in the bottom middle 
panel of Fig. [19] The discrepancies between the best-fit mod- 
els for the blue GCs (shown as thin black lines) and the models 
for the red GCs (thick grey curves) do not permit to prefer any 
specific halo model. 

11.3. Jeans models for a logarithmic potential 

The results are summarised in Table|9] and the model grids so- 
lutions in the (rrj, vo)-plane for the blue GCs (for B - -0.5 and 
0) are shown in Fig. [2T] Again, one notes a strong degeneracy: 
The asymptotic velocity Vo (and hence the total mass) is well 
constrained, while the scale radius ro is not. 

11.4. Joint solutions 

To find a joint solution describing the velocity dispersion pro- 
files of the three tracer populations, we combine the x 2 values of 



the corresponding models and obtain the solution by finding the 
minimum in the co-added^ 2 maps (cf. Sec tl0.6l ). 

1 1 .4.1 . Models for the blue GCs and the stellar velocity 
dispersion profile 

Since the best agreement between models and data can be 
achieved for the blue GCs and the stellar velocity dispersion pro- 
file (see Fig.[T9land the ^-values given in Table[9j, we will first 
combine these two tracer populations to obtain a joint model. 
For the blue GCs, the anisotropy parameters B takes the values 
-0.5, 0, +0.5, while the stellar models are isotropic. The parame- 
ters for corresponding joint models (labelled S.B.tan, S.B.iso and 
S.B.md) are given in Table[9] The best-fit joint (isotropic) mod- 
els are shown in the right panel of Fig.l20l(lower sub-panel). The 
agreement between data and model is best for the two cored halo 
parametrisations: The velocity dispersions of both the stars and 
the blue GCs are very well reproduced by a Burkert halo with 
po = 4.89 x lO~ 2 M pc~ 3 ,ro = 10 kpc or a spherical logarithmic 
potential with ro = 8 kpc and vo = 237 kms -1 . The best-fit joint 
NFW halo, on the other hand, has a very large scale radius and 
over-estimates the velocity dispersion of the blue GCs in the last 
bin (although model and data still agree within the uncertainties). 

12. Discussion 

12. 1. Comparison to the analysis by 
\Chakrabarty & Ravchaud hur\\ 

IChakrabartv & Ravchaudhurvl (I2008I) used the GC kinematic 
database presented in Paper I to study the dark matter con- 
tent of NG C 4636 using the non-para metric inverse algorithm 
CHASSIS (IChakrabartv & Sahall200Tl) . Their main finding was 
that the dark halo required to explain the GC kinematics is very 
concentrat ed. While a high concentration parameter c wlI > 9 as 
derived by IChakrabartv & Ravchaudhurvl is consistent with our 
isotropic Jeans models for the blue GCs (which allow for a wide 
range of concentration parameters), their estimate for the total 
mass exceeds ours: The circular velocity curve shown in their 
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Fig. 19. Observed and modelled GC velocity dispersion profiles. Top row: Models for the blue GCs (sample BlueFinal). From 
left to right, the panels show the best-fit models for an NFW halo, Burkert halo and the logarithmic potential. The solid lines are 
the isotropic models, dashed and short-dashed lines are the tangential (J3 = -0.5) and radial (J3 = +0.5) models, respectively. The 
dash-dotted line is the (isotropic) model without dark matter. The thin vertical line at 1 3 kpc indicates the radial range inside which 
blue and red GCs cannot be distinguished. The data points used in the modelling are shown as filled squares (see also Table 7). The 
model parameters are listed in Table [9] Middle row: The same for the red GCs (RedFinal). Bottom row: Circular velocity curves 
for the best-fit models. Again, from left to right, the results for the NFW halo, Burkert halo and the logarithmic potential are shown. 
The line styles are the same as in the upper graphs, with thin black lines for the blue GCs while the respective models for the red 
GCs are shown as thick grey lines. 



Fig. 6 (left panel) rises to about 450 km s at ~ lOkpc and then 
declines, reaching a value of ~ 370kms _1 at 40 kpc. Our mass 
models (shown in Fig. |T9l > , however, translate to significantly 
lower values of v c with maximal values around 360kms~ 1 (at 
R a 10 kpc) and 300 ^ v c < 340kms~ 1 at 40 kpc. Recently, 



their work has been complemented by an X-ray study which we 
discuss in the following section. 
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Fig. 20. Modelling the stellar velocity dispersion profile. Left: Jeans models for the stars (Bender et al. 1994 data) for an NFW-type 
dark halo and B — 0. The parameters are shown in the (M v ; r , c v ; r )-plane. The thick solid, dashed and thin solid lines indicate the 
68, 90, and 99 per cent confidence limits. The cross indicates the location of the minimum x 2 value. The model parameters are 
listed in T abled The long-dashed (dash-dotted) lines show the median (68 per cent values) for simulated NFW halos as found by 
Bullock et al. (2001). Right: Modelled velocity dispersion profiles for the stars (Bender et al. 1994 data, shown as unfilled squares). 
The upper sub-panel shows the best-fit isotropic models for the stars. The thin solid line shows the model for an NFW-type dark 
halo, and the thin dashed line is the model for the logarithmic potential. The Burkert halo is shown as thick grey line. The thick 
dot-dashed line is the model without dark matter. The lower sub-panel shows the joint models for the blue GCs and the stellar 
velocity dispersion profile (models S.B.iso). The line-styles are the same as in the upper panel. The black squares show the velocity 
dispersion profile for the blue GCs (sample BlueFinal). The halo parameters are given in Table|9] 
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Fig. 21. NGC4636 Jeans models the blue GCs (final sample) where the dark matter component is represented by a logarithmic 
potential Eq. IT6Y Left: Models for a mild tangential bias with B - -0.5. Right: Isotropic models (B - 0). In both panels, the 
best-fit solution is marked by a cross, and the thick solid, dashed and thin solid lines show the 68, 90, and 99 per cent contour levels, 
respectively. The parameters of the halos shown here are also listed in Table|9] 
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Table 9. NGC 4636 Jeans modelling best-fit NFW profiles, Burkert halos and logarithmic potentials 





NFW dark halo 


Burkert halo 


Log. potential 


ID 


Sample 


B 


r s 




Afvir 






X 1 


'"0 


ft 


X' 


*b 


»'o 


X' 








[kpc] 


[M s pc- 3 ] 


[1O' 2 M ] 


[kpc] 






[kpc] 


[M oP c- 3 ] 




[kpc] 


[kms -1 ] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ID 


(12) 


(13) 


(14) 


(15) 


B.tan 


Blue final 


-0.5 


7 


1.01 x 10-' 


1.2 


271 


38.7 


< 0.1 


6 


1.33 X 10-' 


<0.1 


1* 


219 


0.15 


B.iso 


Blue final 





18 


1.47 X 10-- 


2.2 


333 


18.5 


< 0.1 


11 


4.04 X 10- 2 


< 0.1 


4 


234 


< 0.10 


B.rad 


Blue final 


+0.5 


59 


2.07 X 10-- 1 


7.2 


497 


8.4 


0.1 


22 


1.39 X 10~ 2 


< 0.1 


21 


292 


< 0.10 


R.tcm 


Red final 


-0.5 


100* 


5.71 X 10- 4 


6.7 


485 


4.9 


8.9 


100* 


1.94 X 10" 3 


5.4 


100* 


528 


5.2 


R.iso 


Red final 





100* 


6.59 X 10" 4 


8.1 


516 


5.2 


10.9 


100* 


2.13 x 10" 3 


6.5 


100* 


560 


6.2 


R.rad 


Red final 


+0.5 


100* 


8.78 X 10- 4 


11.8 


586 


5.9 


14.4 


100* 


2.81 X 10" 3 


8.4 


100* 


640 


8.0 


S. iso 


Stars (BSG94) 





100* 


9.219 X 10- 4 


12.6 


598 


6.0 


2.2 


100* 


1.90 X 10~- 


1.9 


100* 


1725 


1.9 


Joint solutions for blue GCs and Stars: 


S.B.tcm 


Blue GCs (fi = 


-0.5), Stars (fi = 0) 


100* 


9.66 X 10-" 


13.4 


610 


6.1 


2.2 


8 


6.95 X 10- 2 


2.1 


7 


222 


2.3 


S.B.iso 


Blue GCs (fi = 


0), Stars (B = 0) 


100* 


8.78 X 10" 4 


11.8 


586 


5.9 


2.2 


10 


4.89 x 10" 2 


2.1 


8 


237 


2.2 


S.B.rad 


Blue GCs (fi = 


+0.5), Stars (fi = 0) 


100* 


9.22 x 10~ 4 


12.6 


599 


6.0 


2.3 


13 


3.71 x 10" 2 


2.4 


10 


270 


2.4 



Note 9. In all models, the 7?-band stellar mass-to light ratio is T* ^ = 5.8. Col. (1) labels the models. The second column specifies 
the dispersion profile to which the Jeans models are fit. Col. (3) gives the anisotropy parameter fi. Cols. (4)-(8) lis t the NFW (Eq.fTSl 
and virial parameters (see text for details). Column (9) gives the^- 2 value. The parameters for the lBurkertl(l 1995b halos ([151 , i.e. ro, 
go and the^- 2 value of the best fit model are given in Cols (10-12). Cols. (13-15) are the values for the logarithmic potential (Eq.lTSl. 
Asterisks indicate that the corresponding value is located at the edge of the model grid. 



12.2. Comparison to the analysis bv Uohnson et al.\ 

In their rece n t wor k on the X-ray halo of NGC 4636 
Uohnson et al] d2009l) use a very detailed analysis of 
deep (80 ks) archival Chandra data to derive a mass pro- 
file which they compare to the dynamical modelling by 
IChakrabartv & Ravchaudhurvl (l2008l) . Again, the concentration 
parameters derived for the NFW dark halo models are high, 
with values between 18 and 20. A key finding of their analysis 
is that the derived mass profile depends strongly on whether 
the metal abundance gradient of the X-ray halo is taken into 
account. While the overall shape of the mass profile remains 
the same, the inclusion of the abundance gradient reduces the 
mass at all radii by a factor of about 1.6 (see their Fig. 4). 
Moreover, both models show the same behaviour for large 
radii where th e enclosed mass rises as r 1 2 , a fe ature that was 
also found by lLoewenstein & Mushotzkvl (120031) . To compare 
our dispersion mea surements to the NFW profiles derived by 
iJohnson et alj d2009l) . we proceed as follows: We calculate the 
velocity dispersion profiles expected for the blue GCs for the 
isotropic case (J3 = 0), adopting the NFW parameters given in 
the ir Sect. 4.2. 

IJohnson et alJ parametrise their NFW halos in terms of con- 
centration c and the scale radius r s . TablefTUl lists their values 
together with the corresponding density p$ and the virial param- 
eterU 

For consistency, we adopt for these calculations a distanc e of 
16Mpc, i.e. th e value used by Chakrab artv & Ravchaudhurvl and 
IJohnson et all Using the NFW parameters given in TablefTUl we 
compute the expected velocity dispersion profiles. 

These models for the blue GCs are comp ared to the ob ser- 
vations in the upper panel of Fig. [22] Since IJohnson et all as- 
sumed a very low stellar mass-to light ratio, the difference be- 
tween models Jl and J2, (i.e. the X-ray mass estimate without 
abundance gradient) before and after the subtraction of the stel- 
lar component is small. The corresponding velocity dispersions 
lie well above the data points. 

A much better agreement between X-ray and GC based mass 
estimates is achieved when the metal abundance gradient of the 

6 Note that Johnson et al. ( 2009) use a different definition of the virial 
parameters, i.e /?2uo = c - r„ where i? 2 oo is the radius within which the 
mean density equals 200 times the critical density of the Universe. 



Table 10. Parameters of the NFW halos derived bv lJohnson et all 
(2009) 



Model 


T s 


c 


^200 




M 2m 




[kpc] 




[kpc] 


[Mope" 3 ] 


[1O 12 M ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Jl 


21.8 + 0.9 


20.1+0.8 


438 + 25 


0.0359 + 0.0036 


9.8! 2 ;* 


J2 


24.6 + 0.9 


18.0 + 0.6 


443 + 22 


0.0270 ± 0.0023 


10 -°-L9 


J3 


26.1 ± 1.0 


14.4 + 0.4 


376 + 18 


0.0154 + 0.0011 


6 2 +u 



Note 10. T he first Col, labels the models in order of their appearance in 
Sect. 4.2 of Joh nson et all {2009), where Jl: X-ray data, total mass; J2: 
X-ray data, stars subtracted; J3: X-ray profile including metal abun- 
dance gradient, stars subtracted. Columns 2 and 3 are the parameters 
quoted by these authors, Col. 4 gives R 2 m in units of kpc. The corre- 
sponding values for the density g s and M200, i.e. the enclosed mass at 
^200 are given in Cols. 5 and 6. 



X-ray gas is taken into account. Model J3 agrees, within the 
uncertainties, with the GC data out to about 30 kpc. For the 
abundance gra dient corrected mass profile shown in Fig. 4 of 
John son et all one obtains a very similar velocity dispersion pro- 
fiidZl The r rise for large radii, however, leads to an almost 
constant velocity dispersion profile for R ^ 40 kpc. For refer- 
ence, we also plot, in Fig. [22] the best-fit isotropic NFW model 
derived for the blue GCs assuming a distance of 16 Mpc and the 
adjusted Y* « = 6.4. 

12.3. Are all GCs bound to NGC 4636? 

Objects with velocities in excess of the escape velocity are prob- 
able interlopers. Due to the logarithmic divergence of the NFW 
potential, the escape velocity is not defined. But, in any spherical 
potential bound particles travel on planar orbits, and energy and 
angular momentum conservation are used to derive the following 
expression: 

2 _ 2r fl 2 (cD(r p )-4>(r fl )) 



7 For this calculat ion we used a piecewise fit to the data in Fig. 4 of 
IJohnson et all J2009h 
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where v p is the pericentric velocity, r p and r a are the pericentre 
and apocentre distances, respectively. 
The gravitational potential (D(r) given by 



O(r) = -4tt 



— I £>(s)s 2 ds + I g(s)sds 
r Jo Jr 



(18) 



where g = p stars + £>dm is numerically integrated using the NGC 
4636 stellar mass profile and the NFW halo dark matter density 
profile obtained from the blue GCs (model B.iso). 

Objects outside a given curve have apocentric distances 
larger than the corresponding value of r a . The set of 
curves shown in Fig. [23] is obtained from Eq. [T7] by fixing 
r a e {40, 60, 100, 150, 200, 300}kpc) 

For the two blue GCs (objects 3.1:69 and 3.2:65) with good 
velocity measurements (Vheiio = 1428 + 37 and 1441 +28 km s , 
respectively) at a galactocentric distance of as 34 kpc, we find 
apocentric distances of more than 150 kpc. Given that these con- 
ditions are extreme, an unknown population of GCs with large 
apogalactic distances may be present in the bulk of velocities. 
However, the question whether these GCs are bound or unbound, 
cannot be answered. Recall that, even in the Milky Way sys- 
tem, some GCs have Galactocentric distances of more than 100 
kpc. How do these objects compare to the GCs with surpris- 
ingly high relative v elocities in the NGC 1399 GCS identified by 
Richt ler et al.l (120041) ? These authors show, in their Fig. 20, that 
the objects in their sample of (about a dozen) GCs with heliocen- 
tric velocities below 800 km s (which corresponds to velocities 
of at least 640 km s _1 with respect to NGC 1399) have apogalac- 
tic distances between 100 to 200 kpc (with one GC even featur- 
ing ra A 400 kpc). However, as shown in lSchuberth et al.l (12008). 
blue, velocity-confirmed NGC 1399 GCs are still found at dis- 
tances of about 200 kpc. Thus, while the large number of GCs 
with apogalactic distances of r a m 200 kpc might be surprising, 
it is well possible that these objects bel ong to the very extend ed 
NGC 1399 GCS. Another scenario (cf. [Schuberth et al.ll2010l) is 
that these metal-poor GCs were stripped from infalling galaxies 
which is not unlikely for a galaxy such as NGC 1399 which is 
the central galaxy in a relatively dense cluster. In the NGC 1 399 
GCS, the most extreme combination of radial di stance and ve- 
locity is that of gc381.7 (from the catalogue of Bergondet al. 
2007) which would have an apogalactic distance of the order 0.5 
to 1 Mpc (Schuberth et al. 2008), i.e. of at least twice the Fornax 
cluster core radius. 

In the case of NGC 4636, the two blue GCs with extreme 
velocities are remarkable. Although the estimated apogalactic 
distances are smaller than those of the extreme GCs near NGC 
1399, one has to take into account that the NGC 4636 GCS ap- 
pears to be truncated, and that almost no GCs are found beyond 
as 60 kpc. Moreover, NGC 4636 is relatively i solated and doe s 
not show any signs of recent major mergers (iTal et al.l 12009). 
This would make these objects candidates for a population of 
'vagrant' GCs belonging to the Virgo cluster rather than being 
genuine members of the NGC 4636 GCS. 

12.4. Stellar mass-to-light ratio 

Our value of the stellar M/L^-value of 5.8, which we need under 
isotropy to model the central velocity dispersion data, is not di- 
rectly supported by theoretical single stellar population synthe- 
sis models. Depending on the adopted stellar initial mass func- 
tion (IMF), mod els predict values around 3.5-4 for metal-rich 
old populations {Bruzual & Chariot! I2003I: iThomas et al.l I2003 : 
Mara ston et al.l l2003b iPercival et al.ll2009l) and can reach up to 



4.7 for super-solar abundance even for IMFs with a gentle slope 
in the mass-poor domain like the Kroupa (2001) IMF. 

Of course, one could surely find an appropriate radial bias 
which at a given mass enhances the projected velocity dispersion 
in the central regions and thus would permit a lower M/L-value. 
This sort of fine-tuning is somewhat artificial and moreover is 
not supported by observational evidenc e. NGC 4636 exhibits in 
the analysis of lKronawitter et al. (2000) a tangential bias, though 
at larger radii, and is isotropic in its inner region. 

A good reference to stellar Mj L -v alues of the inner region s 
of elliptical galaxies is the study of Cappella ri et al.l ((2006). 
These authors used the SAURON integral-field spectrograph to 
compare the /-band dynamical M/L with the (M/L pop ) obtained 
from stellar population models for a sample of 25 early-type 
(E/S0) galaxies. 

For the 24 sample galaxies which lie in the /-band magnitude 
range -20 5 M/ $ -24 they find a correlation between the dy- 
namical M/L and the galaxy luminosity, in the sense that M/L 
weakly increases with luminosity as L 7 ° 32 (see. their Fig. 9, and 
Eq. 9 for the fit). 

To address the question whether these observed M/Ld yn vari- 
ations are due to a change in the s tellar populations or to differ- 
ences in the dark matter fraction, Cappella ri et al.l plot in their 
Fig. 17 the dynamical M/Ldyn as a function of M/L pop . Both 
quantities are correlated but while M/Ldyn ^ M/L pop for all 
sample galaxiefl the data points clearly lie off the one-to- 
one relation. The authors consider the lower luminosity fast ro- 
tators and the high luminosity slowly rotating galaxies sepa- 
rately (NGC 4636 would belong to the latter group). For old (age 
> 7 Gyr) galaxies they find that the dark matter fraction within 
one effective radius increases from zero to about 30 per cent as 
the dynamical mass to light ratio increases from 3 to 6. At a given 
(M/L pop ), the massive slow rotators have higher dynamical M/L 
values than the less massive fast rotators. How does NGC 4636 
fit into this picture? 

For NGC 4636 the iMaraston etaD (120031) SSP model pre- 
dicts an M/L p0 pj = 3.27 (for solar metallicity and an age of 
13 Gyr). Converting this to the B and R bands, one obtains 
M/L p0PiB = 6.6 and M/L p0PjS = 3.9, respectively. Gerhard et 
al. quote an even lower value M/L p0Pv B = 5.9. 

What woul d we expect from the relations given by 
Cappella ri et alJ ? 

From Table 1 , one obtains M jjq% = -23.3, and (for M?, Q = 

4.08 Langl d!999l) ) from Eq. 9 in Cappellari et all we thus would 
expect a dynamical M/L/ = 4.7, corresponding to M/Ld yn R = 
5.7. 

However, the models of Cappellari et al 1 (120061) have.bv def- 
inition, a radially constant M/L, while our M/L depends on ra- 
dius and reaches M/ Ld vn „R ~ 8 at the effective radius. 

Cap pellari et al.l speculate on the possibility that the differ- 
ence between dynamical and population M/L is due to a higher 
dark matter content of more luminous galaxies, but the general 
question is whether it is appropriate to apply SSP models to 
composite stellar systems. Let us consider u> Centauri, probably 
the dynamically best investigated stellar syste m, which is unre- 
laxed and composed of different populations. Ivan de Ven et al.l 
(2006) quote a V-band M/L of 2.5 ±0.1. The metallicity dis- 
tribution of stars in u Cen has a maximum at [Fe/H1»-1.7 
with a broad tail towards higher metallicities (e.g. iHilker et all 



8 M /L DOD was estimated using the Vazdekis et al. (1996); Vazdekis 
( 1999) stellar population models assuming a lKroupal d200lh IMF, and 
the dynamical M/L was obtained from Schwarzschild modelling of the 
SAURON data. 
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2004; Calami da et all 12009). The more metal-rich populations 
are probably also young er by a few Gyr. F rom the population 
synthesis market, we cite lPercival et al.l (12009b who quote 2.3 as 
the value for a population with [Fe/H]= —1.7 dex and an age of 
13.5 Gyr, and 2.0 for a population with [Fe/H]= -1 .3 dex and an 
age of 10 Gyr, adopting a Kroupa IMF. Without aiming at pre- 
cision, the composite 'population' M/L will probably not reach 
the dynamical value of 2.5, unless there are old metal-rich pop- 
ulations, for which there is no evidence, so a> Centauri is at least 
a mild example without dark matter, where the dynamical mass 
is larger than the population mass. 

However, an elliptical galaxy is a composite system with a 
long and complicated star formation history. If the IMF in a lo- 
cal star formation event is universal, there is no guarantee that 
the final mass function in a galaxy bears the same universality. 
Star formation occurs in star clusters and a galaxy's field popula- 
tion is composed of dissolved star clusters. If the mass spectrum 
of star clusters is a power-law like rtT 1 then the dissolved pop- 
ulation is the result of adding up many low-mass clusters, but 
fewer high-mass clusters, where the full ste llar mass spectrum 
can be expected. Weidner & Kroupal J2006) showed that, if the 
maximum stellar mass within a cluster depends on the clusters' 
mass, the resulting stellar mass function can be even steeper than 
a Salpeter mass function. There are no simulations of the final 
M/L of an elliptical galaxy available, but since a Salpeter-like 
mass function increas es the M/L by factor of roughly 1.4 (e.g. 
Cap pellari et al.l2006b . it is plausible that there is not a strict uni- 
versality of stellar mass functions among galaxies, but that the 
stellar mass function of an old elliptical galaxy may depend on 
the history of its assembly. In conclusion, a stellar M/L^-ratio of 
5.8 might well represent the stellar population. 

Another consideration may be worthwhile: 

If we require the M/L-values to agree with the SSP predic- 
tions, we need a M/Lr = 4 or smaller, lets say, 3.7. The dark 
halo, represented by a logarithmic halo, would assume parame- 
ters like ro = 1 kpc and vo = 250 km/s. The central density of dark 
matter then is 3.5 M /pc 3 under isotropy, and equality of stellar 
mass and dark mass is reached already at a radius of about 3.5 
kpc. The central projected velocity dispersion is 170 km/s for the 
stellar mass alone and 192 km/s for the total mass. If that would 
be typical for elliptical galaxies (of which there is no evidence), 
scaling relations like the fundamental plane would dynamically 
be dominated by dark matter and the 'conspiracy' between dark 
and luminous matter would reach a level even more difficult to 
understand than it is now. 

Finding such a high central dark matter densi ty prompts us 
to con sider an older argument brought forward by Gerhar d et alj 
(2001): the dark halos of elliptical galaxies in their sample turned 
out to exhibit a central density which is higher by a factor of 
at least 25 than those of spiral galaxies of similar luminosity, 
and also that the phase space densities are higher. Since in col- 
lisionless mer ging ev ents phase space densities cannot grow, 
iGerhard et alj argued that it is unlikely that dark halos of ellipti- 
cals formed by the merging of dark halos of present-day spirals. 

Gerhard et al.'s expression for the phase space density reads 
//, = 2 3 l 2 ph/v 3 h , being p/, the central density and v? the character- 
istic halo velocity of a logarithmic potential. 

The above hypothetical dark matter density i s about a fac- 
tor 10 00 higher than that of spirals (see Fig. 18 of iGerhard et alj 
1200 ll) and the phase space density of the corresponding halo is 
6.7 xl0~ 7 in units solar masses, pc, km/s, much higher than those 
of spirals. We conclude that with our example low M/L, it might 
not be possible to reach these densities and phase space densities 



by collionless accretion of spiral-like halos, and one has to resort 
to dark halos resembling those of dwarf spheroidals. 

12.5. MOND related issues 

The question whether elliptical galaxies fulfill the predictions 
of Modified Newtonian Dynamics (MOND, see e.g. Milgrom 
2009, Sanders & McGaugh 2002) obviously is a fundamental 
one. Ellipticals have so far been less in the focus of MOND 
than disk galaxies. The most compelling case of an appar- 
ently MONDian elliptical galaxy is the E4 galaxy NGC2974 
( Weii mans et al.ll2008h where the extended Hi disk permits a se- 
cure determination of the circular velocity which is constant out 
to 20 kpc. 

In Paper I, we noted already that NGC 4636 seems to be con- 
sistent with being MONDian. Here, we plot, in the lower sub- 
panel of Figf24] the MOND circular velocity curve (shown as 
thick solid line) for the stellar mass profile of NGC 4636 (for 
M/L +> r=5.8). The MOND circular velocity curve is obtained 
from the Newtonian one via the following equation: 



V 1 



v 4 



+ V; 



irc,N 



a ■ r 



(19) 



where V c ; 1c n is the N ewtonian velocity. F or ao we adopt the 
value recommended by Fam aev et al.l J2007): 1.35 x 10 -8 cms~ 2 . 
Within the central « 40 kpc, i.e. the radial range, for which we 
have data, the MOND circular velocity curve agrees fairly well 
with the best-fit joint model for stars and the blue GCs (models 
S.B.iso, Burkert halo, shown as thin dashed line). 

Below we put NGC 4636 into the context of the more recent 
literature. 



12.5.1. Comparison to the /zoo-relation bv lDonato et alj 

Recently Donato et alJ d2009|) . in their extension of the work by 
iKormendv & Freem an (2004), confirmed that the central surface 
density of galaxy dark matter halos is nearly constant and al- 
most independent of galaxy luminosit y. MOND-related as pects 
of this finding have been disc ussed by [Gentile et al.l (I2009I) and 
iMilgroml (I2009I) . Donato et all assume that the dark matter halos 
of the galaxies are described bv lBurkerfl(ll995l) halos (cf. Eq.[T4l 

and[T5)- 

For the DM surface density //od = rogo, iDonato et alJ find 
the following relation to hold for galaxies in the magnitude range 
-8 > M B > -22: 



log 



Hod 
M e pc- 



= 2.15 ±0.2 



(20) 



How does NGC 4636 fit into this picture? From the values 
given in Table[9] it appears that the dark matter density of our 
Burkert halo models (for the blue GCs and the stars) (2.50 ^ 
log(/i()D / M e pc~ 2 ) S5 2.90) is too high with respect to the above 
value. 

Going back to Fig. [18] however, one sees that for ro ^ 20 kpc 
the 68 per cent CL contour (thick solid line) of our isotropic 
models for the blue GCs lies within the range of values from 
Eq. [20] (shown as dot-dashed lines). In the upper panel of 
Fig. [24] we show, as an example, isotropic Jeans models for 
the blue GCs and the stellar component for ro = 20 kpc, p p = 
1.115 x 10~ 2 , logyL/o.o = 2.35 (solid lines) , i.e. a Burkertfhalo 
whic h is consistent (wit hin the uncertainties) with the relation 
from [Donato et al.ll2009l For reference, the best-fit joint model 
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is shown with dashed lines. The lower panel of Fig. [24] shows 
the corresponding circular velocity curves. Regarding the stars, 
we find^- 2 values of;f 2 oint = 2.0 an d^ onato = 2.9 for the joint 
model and the model which is consistent with Eq. [20] respec- 
tively. For the blue GCs, the respective values are xf omt < 0.1 

andy 2 , =0.7. 

^ donato 

We conclude that the GC and stellar dynamics of NGC 4636 
do not contradict the constant dark matter density relati on, but 
we caution that the halos of ellipticals in the lDonato et al J sample 
are constrained only by weak lensing shears, which do not probe 
the inner regions. 

iGerhard et al.L for their sample of ellipticals, implicitly 
found (multiplying their equations (6) and (8)) po ■ ro = 
635-Lj"j/jo.65. which for all practical purposes is constant (they 
used logarithmic halos instead of Burkert halos). This fits better 
to our value and the question, whether disk galaxies and ellipti- 
cals really show the same surface density of dark matter, remains 
open. 

12.5.2. The baryonic Tully-Fisher relation and NGC 4636 

A severe irritation to the ACDM paradigm on galactic scales is 
the existence of a baryonic Tully-Fisher relation (BTFR) among 
spiral galaxies with an astonishingly small scatte r, covering five 
orders of magnitude in mass, which reads (lMcGaughl2005l) : 

M bar = 50 ■ v* at , (21) 

where vg at is the circular velocity (in units of kms~') at a radius 
where the rotation curve becomes flat, and Mb ar is the total mass 
in baryons (in units of M Q ). Such a relation would nat urally re- 
sult from Modified Newton ian Dy namics (MOND, e.g. lMilgroml 
ll983LISanders & McGaughll2002l) . In the deep MOND regime: 

4 

VZ ■ an 

M = - — - , (22) 
G ' 

M being the total mass, G the gravitational constant, and an the 
MOND constant, w hich has the value 1.35^-^ x 10~ 8 cms~ 2 
(Famaevetal. 2007). The factor of 50 in the above relation 
(Eq. |2TT i corresponds to a somewhat higher value of an = 
1.5 x 10~ 8 cms~ 2 which still lies within t he uncertainties of 
the value quoted above. The recent work of St ark et al.l (f2009) 
and Trachternach et al. (2009) confirmed this relation, which in 
IStark et al.l formally reads Mb ar = 61 • Vg at 4 agreeing even better 
with the canonical value for an. 

There als o exists a BTFR for elliptical g alaxies 
(IGerhard et al l 120011; iMagorrian & Ballantvnel 1200 It 
Thomas et al. 200 7]), but most el liptical galaxies lie off the 
spiral relation (IGerhard et al.ll200"Tl) . Since the outermost radii in 
these dynamical studies may be still on the declining part of the 
circular velocity curves, it is interesting to put NGC 4636 into 
this picture. Although we cannot strictly distinguish between 
a constant circular velocity curve and a slightly declining one, 
the model with the flattest rotation curve has a circular velocity 
of 300kms _1 and thus we expect a total baryonic mass of 
4 x 10 n M o . With the data from Table 1, we have M R = -22.6, 
thus 5.5 X 10 10 L G and a total baryonic mass of 3.2 x 10 n M Q , 
which would place NGC 4636 a bit below the relation for 
spirals. However, given the uncertainties in adopting distances, 
M/L-ratio and even the absolute solar R-magnitude, we are 
reluctant to assess these value as a clear displacement and repeat 
our conclusion from Paper I that NGC 4636 is consistent with 
the MONDian prediction. 



In any way it is of fundamental interest to investigate more 
elliptical galaxies at large radii. If ellipticals and spirals would 
follow the same BTF-relation in spite of very different formation 
histories, the challenge for the Cold Dark Matter paradigm of 
galaxy formation would be considerable. 



13. Conclusions 

We revisit the dynamics of the globular clusters system of NGC 
4636 on the basis of 2 89 new globular cluste r (GC) velocities. 
Including the data from l Schuberth et aLl (I2006I Paper I), our total 
sample now consists of 460 GC velocities, one of the largest 
sample obtained until now for a non-central elliptical galaxy. In 
addition we present new kinematical s t ellar d ata extending in 
radius the analysis bv lKronawitter et aLl (l2000l) . 

We model the total mass profile by the sum of the stellar 
mass plus a dark halo, for which we adopt different analytical 
forms. With our distance of 17.5 Mpc, we need a stellar M/Lr- 
value of 5.8 to satisfactorily reproduce the projected stellar ve- 
locity dispersion near the center under isotropy. This value is 
higher than values predicted from canonical population synthe- 
sis of an old, metal-rich popu lation, resembling the results from 
the SAURON collaboration (ICappellari et al.ll2006b . We argue, 
however that the actual stellar mass function of an elliptical 
galaxy might be somewhat steeper than IMFs of local young star 
clusters. 

We perform spherical Jeans-analyses independently for the 
red and the blue cluster population and fit dark matter profile 
parameters for NFW-profiles, for Burkert profiles and for log - 
arithmic halos for different anisotropies. The fits using the red 
cluster populations are consistently worse than using the blue 
populations, a finding, which differs from our previous study of 
the central cluster galaxy NGC 1399 and for which we do not 
have a good explanation. Our recommended joint logarithmic 
halo which uses the stellar light and the blue GCs, has the pa- 
rameters rn — 8 kpc and vo = 237 km/s. 

The higher moments of the velocity distributions of the blue 
and red GC subpopulations are not really stable against the sam- 
ple selections and thus do not permit to seriously constrain pos- 
sible GC orbit anisotropies. However, they are consistent with 
the GC orbits being isotropic to a good approximation. 

We compare ou r results with the mass profile, derived from 
X-ray analysis, of iJohnson et al.l (120091) . When the element 
abundance gradient of the X-ray gas is not taken into account, 
the X-ray mass profile exceeds the GC mass profile by a sig- 
nificant factor. If the abundance gradient is accounted for, the 
agreement is good out to 30 kpc, but the X-ray mass profile still 
exceeds our mass profile beyond this radius. The might be a gen- 
eral problem of X-ray analyses, if strong abundance gradients 
are present, for example in galaxy clusters. 

NGC 4636 almost falls onto the baryonic Tully-Fisher rela- 
tion for spirals and behaves more or less MONDian, as already 
noted in Paper! However, when modelled with a logarithmic 
halo, its halo surface density resembles that of other elliptical 
galaxies. 
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Appendix A: The velocity dispersion profile of 
NGC4636 

Appendix B: The stellar mass profile of NGC4636 

Here we give the piecewise approximation of the stellar mass profile of 
NGC4636 (see Sect. 1 10.21 used in our modelling. In the following expressions, 
x is in units of parsec, D is the distance of NGC4636 in Mpc (assumed to be 
D=17.5 in our modelling), and r (= 5.8) is the R-band stellar mass-to-light 
ratio. 

forO < x < 3.0 ■ D: 



/ D\ 2 

ft = Mb) - 5J4 



x 10° 



for 3.0 ■ D < x < 29.184 ■ D: 



(B.l) 



(B.2) 



whereao = 1.322xl0 7 ,fli = -1.118x 10 7 , a 2 = 3. 158x10 s , a 3 = -9.127xl0 4 , 
anda 4 = 1.164 X 10 3 . 

for 29.184 ■ D< x < 104.49 ■ D: 



(B.3) 



where b\ = 9.616 x 10 6 , b 2 = 1.034 x 10 6 , 63 = -9.133 x 10 3 , and b 4 = 25.78. 
for 104.49 • D < x < 2597.9 ■ D: 



P3(x) = T 



(B.4) 



where c = 6.101 X 10 s , ci = 4.31 X 10 7 , c 2 - 1.507 X 10 4 , and c 3 = 1.996 
for x > 2597.9 ■ D 



A(x) = d\ ■ arctan 



d Q ) + dl ( do ) +dl ( d 



(B.5) 



where 

d = 1.194 X 10 3 ,rfi =4.221 X 10 l0 ,d2 = -0.0668, and d3 = 0.003774. 

Appendix C: Globular cluster velocities 
Appendix D: Foreground star velocities 
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Table A.l. NGC4636 velocity dispersion profile from FORS 2 
data 
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Note 11. Column (1) gives the slit number on the FORS 2 Mask 1_1 
from Paper I. The galactocentric distance R (in units of arcseconds) is 
given in the second column. Positive and negative values of R refer to 
positions to the north and south of the centre of NGC 4636 respectively. 
Column (3) lists the velocity dispersion (in km s~'). 
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Fig. 12. Velocity dispersion profiles. Upper panel: Circles show 
the dispersion values for the blue GCs (for fixed radial bins 
cf. Table[6} prior to quality selection and outlier removal. Dots 
are the values for the final blue sample. The dashed horizon- 
tal 'error bars' indicate the radial range of a given bin; the bins 
used here are the ones indicated by the dotted lines in the lower 
panels of Fig [8] The dispersion values from Paper I are shown 
as diamonds. In both panels, the vertical line at 2.'5 indicates 
the first bin which covers the radius inside which blue and red 
GCs cannot be distinguished. Lower panel: The same for the 
red GCs. 
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Fig. 13. The stellar velocity dispersion profile of NGC4636. 
The diamond is the central dispersion measured by 
(120021) . Large dots show the values from 
squares are the values derived from our 
Unfilled (filled) squares 
The vertical 



994), 



Proctor & Sansor 
Bender et all (11994 
FORS 2 spectra (see text for details), 
indicate slits North (South) of the galaxy centre, 
dashed lines indicate the region we use to model the galaxy's 
velocity dispersion profile (cf. Sect. II 1.1.31 1. The solid green 
curve is a linear fit to the data fo r R > 0.'5. The major axis 
measurements bv lPu& Hani (1201 lh are shown as small (orange) 
dots. For radii < 0.219', where the long slit extends to either 
side of the galaxy, we plot the averaged values; the uncertainties 
were added in quadrature. For comparison, the data points for 
blue and red GCs for constant number bins (final samples, 
values from Table|7]i are shown as triangles. The vertical dotted 
line at R — 2'.5 marks the radius inside which blue and red GCs 
cannot be separated. 
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Fig. 14. Deprojecting the surface brightness profile of 
NGC4636. Upper sub-panel: /?-band surface brightness 
profile. The data points are from D05+ (their Table A.4), and 
the dashed line is the fit given by D+05 (Eq. [7j. The solid line 
is the three-component fit given in Eq. [8] and the thin solid 
lines show the individual components. Bottom sub-panel: 
Luminosity density profiles in units of L pc~ 3 . The dashed line 
was obtained by numerically integrating Eq. Q The solid line is 
the analytical deprojection of Eq. [8j again, the components are 
shown as thin solid lines. The vertical dotted lines indicate the 
radial range of the data points shown in the upper panel. 
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Fig. 15. NGC 4636 stellar mass profile. The thin solid curve is 
Eq. [10] for Y* « = 5.8 (obtained through numerical integra- 
tion). The thick dotted curve is the piecewise approximation 
used in our modelling. The radial range of the respective pieces 
is indicated by the bars, and the functions and coefficients are 
given in AppendixlBl The vertical dashed lines indicate the ra- 
dial range of the photometric data by D+05. The horizontal line 
at 6.15 x 1O U M shows the asymptotic value, i.e. M(r = oo). 
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fdark [kpc] 

Fig. 18. Jeans m odels for the blue GCs (sample BlueFinal, ft = 
0) for a iBurkertl dark matter halo. The cross marks the best- 
fit parameters (cf. Table|9}- The thick solid, dashed and thin 
solid lines indicate the 68, 90 and 99 per cent confiden c e lim- 
its. The thick dot-dashed line indicates the lDonato et al.l (2009) 
central dark matter surface density relation logyuoo = 2.15 ± 
0.2 [log(M Q pc~ 2 )] (the thin dot-dashed lines show the uncertain- 
ties). 




Fig. 16. Number density profile of the blue and red GCs. The 
data from D+05 (their Table A3) are shown as crosses and dots 
for the blue and red GCs, respectively. The curves are the fits 
(cf. Eq.Qj] Table© to the data, and the radial range 2 .'5 < R < 
13.'5 is indicated by the arrow and the vertical lines. 
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Fig. 22. Comparison to the NFW halos derived bv lJohnson et al.l 
(2009) with the parameters listed in Table[10] Upper panel: 
Velocity dispersion profiles. The data points show our final blue 
GC sample (same as in Fig. [19] left panel but for a distance of 
16Mpc). The thick solid line shows model Jl; model J2 is shown 
as short-dashed line, the dashed regions show the corresponding 
uncertainties. The long-dashed line is model J3, and the grey 
area shows the uncertainties. The dash-dotted line (labelled X^) 
corresponds to the mass pro file (incorpora ting the abundance 
gradient) shown in Fig. 4 of iJohnson et al.l The thin solid line 
is the best-fit model for the GCs (r s = 20 kpc, p s = 0.012, 
T*,« = 6.4). Lower panel: Circular velocity curves. The line 
styles are the same as in the upper panel. 
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Fig. 23. Radial velocity versus galactocentric distance for blue 
GCs with velocity uncertainties Av < 50kms _1 . The curves, 
calculated from Eq.[l7]for the total mass as given in model B.iso 
(Table[9] NFW halo), indicate pericentric velocities for fixed 
apocentric distances. Objects outside a given line have apocen- 
tric distances of at least the value for which the curve was calcu- 
lated. 
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Fig. 24. Comparison to the results from iDonato et al.l (|2009). 
Upper panel: Open squares show the stellar velocity dispersion 
profile (BSG94), filles squares are the blue GCs (BlueFinal). The 
velocity dispersion profile for the joint isotropic model (blue 
GCs and stars, S.B.iso, Burkert halo) is shown as dashed line. 
The solid line shows the dispersion for a Burkert dark matter 
halo ro = 20 kpc, go = 1.115 X 10 ~ 2 (i.e. log^np _ 2. 35) which 
is consistent with the relation by Do nato et al.l (120091) (see text 
for details). Lower panel: Circular velocity curves. Again, the 
thin solid line shows the curve for the halo with log 0D = 2.35, 
and the dashed line is the joint model. For reference, the circu- 
lar velocity curve for the stars alone (Y*s = 5.8) is shown as 
long-dashed line. The MOND circular velocity curve (Eq.[T9T> is 
shown as a thick solid line. The dotted vertical line at 40 kpc in- 
dicates the location of the outermost velocity dispersion bin for 
the blue GCs. 
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Table D.2. Catalogue of globular cluster velocities. Column (1) gives the spectrum identifier (mask:slit). Note that the new data have the prefix 
"f". For the purpose of cross-referencing, the data from Paper I are also included in this table. Columns (2) and (3) are the right ascension 
and declination, respectively (for objects with Washington photometry, the coordinates are those from Dirsch et al. 2005; for objects not in 
the photometric catalogue, the values were estimated from the pre-images). The projected distance from the centre of NGC4636 (in units of 
arcminutes) is given in Col. (4). Column (5) is the position angle, (measured North over East). Columns (6) and (7) are the magnitude and colour 
as given by Dirsch et al. (2005). For the new objects, Column (8) gives the heliocentric velocity measured for a given spectrum, for the data from 
Paper I, the final velocity given there is quoted. Column (9) gives the final velocity used in this study, for objects for which more than one spectrum 
was obtained, the error-weighted mean is given. Column (10) indicates the subsample ("R.25" denotes GCs with projected distances R < 2.'5). 
Column (11) gives the cross-identifications within this list. The entry in Col. (12) is T for the first appearance of an object in this table. All 
subsequent entries for the same object are labelled '0'. 
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RedFinal 
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133.5 
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RedFinal 
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134.5 


21.27 
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Table D.2. continued. 
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,09.10. m 9 
+U2. Jo. JU.2 


J.JO 


1 AS A 
14J.4 


1 Q 81 

ly.oj 


1 77 
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1 AiQ-i- 1 m 
1Ujo± 1U / 


Blue 


f A 1 A O 

iu 1 -4y 


U 


IUZ-Z I 


1 7-zll-n 1 (\S 
1Z.4J.U1.UJ 


+uz. jy.u4. / 
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Table D.3. Catalogue of foreground star velocities. Column (1) gives the spectrum identifier {mask: slit). Note that the new data have the prefix 
"f". For the purpose of cross-referencing, the data from Paper I are also included in this table. Columns (2) and (3) are the right ascension 
and declination, respectively (for objects with Washington photometry, the coordinates are those from Dirsch et al. 2005; for objects not in 
the photometric catalogue, the values were estimated from the pre-images). The projected distance from the centre of NGC4636 (in units of 
arcminutes) is given in Col. (4). Column (5) is the position angle, (measured North over East). Columns (6) and (7) are the magnitude and colour 
as given by Dirsch et al. (2005). For the new objects, Column (8) gives the heliocentric velocity measured for a given spectrum, for the data from 
Paper I, the final velocity given there is quoted. Column (9) gives the final velocity used in this study, for objects for which more than one spectrum 
was obtained, the error-weighted mean is given. Column (10) gives the cross-identifications within this list The entry in Col. (1 1) is '1' for the first 
appearance of an object in this table. All subsequent entries for the same object are labelled '0'. 
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Table D.3. continued. 
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- IZ / . j 


70 7Q 

zu.3y 


7 7A 

Z. 10 


1 1 O-i-C;^ 
-1 lU+JJ 


1 1 n_i_ ^ ^ 
-1 1U+JJ 






7 7-7C 
Z.Z. JJ 


1 7 - /1 7 - 7 n 11 

1Z.4Z.jU. 1 1 


_i_m-7A-77 A 

+UZ. jO.Zj.O 


A O/l 

o.y4 


1 7/1 A 
-134.0 


1 se 

iy.00 


1 77 

1 .3Z 


1 7S_i-/1 ^ 
-130+4J 


1 78-1./1 

-130+4J 






2.2:50 


12:42:16.37 


+02:37:34.8 


9.14 


-113.8 


19.08 


0.93 


-42±40 


-42±40 






2.2:53 


12:42:19.87 


+02:37:48.9 


8.25 


-114.7 


19.80 


4.30 


-48±85 


-48±85 






2.2:55 


12:42:18.69 


+02:38:02.2 


8.43 


-112.5 


18.49 


0.54 


-67±26 


-67±26 






2.2:66 


12:42:29.58 


+02:38:46.8 


5.64 


-116.1 


21.15 


3.08 


-188+65 


-188±65 






2.2:68 


12:42:21.32 


+02:38:54.0 


7.51 


-108.4 


20.17 


1.14 


237±22 


237±22 






2.2:75 


12:42:28.04 


+02:39:23.6 


5.76 


-109.0 


19.97 


3.36 


-31 ±47 


-31 ±47 






2.2:76 


12:42:29.86 


+02:39:27.1 


5.32 


-110.0 


21.90 


1.40 


240±39 


306±6() 


f 10-06 






Table D.3. continued. 



ID 


R.A. 


Dec. 


— p 

"4636 


PA 




C-R 


^spec 




r^rn*;*;-Trlpntif ip?itinn^<;^ 


1 st 




yJLUVyJ) 


/T9000~i 
(^J zuuu ^ 


[arcmin] 




[mag] 


[mag] 


[Km s j 


[~L"m 

[Km S J 






9 7-7Q 

L.L. to 


1 7-/17-70 7C 

1z.4z.zy. jj 


l Q9-7Q-7C 7 

+uz.jy.jj. / 


C 7Q 

j.jy 


1 ne 1 
-1U5. 1 


1 8 A/1 
16.04 


n 77 
U.JJ 


Q7_i_9Q 

-y /±zo 


no j. oo 

-y /±zo 




1 
1 


7 7-SQ 

z.z.sy 


1 7-47-77 QQ 

1z.4z.zj.yy 


0-07-40- 1 7 A 
+UZ. 4U.1 / .0 


A C.4 
0. J4 


QS A 

-yo.o 


7 1 AQ 

zi .oy 


O 7S 
U. / J 


97A-1-7Q 
-ZjO± JO 


77A-1-78 
-ZJO± Jo 




1 

1 


o 7-qc. 
z.z. yj 


1 7-4.7-71 OS 
1Z.4Z.J1.U5 


, H7-40-7A 4 
+UZ.4U. JO. 4 


4 74 
4. /4 


QS O 

-yo.u 


70 7A 
ZU. JO 


9 CO 
Z. jZ 


74-1-47 
J4±4Z 


74-1-49 
J4+4Z 




1 
1 


7-QA 

z.z. yo 


1 7-47-70 87 
lZ.4Z.ZV.oJ 


_i_07-40-40 7 
+UZ.4U.4U. J 


C. 04 
J.U4 


QA S 

-yo.o 


70 A1 
ZU. 01 


7 94 
J.Z4 


1 7-1-44 
-1 /±44 


1 7-1-44 
-1 /±44 




1 
1 


7-Q7 

z.z.y / 


1 7-47-70 01 

1z.4z.zy.u1 


, 07-40-44 1 
+UZ.4U.44. 1 


K 94 
J.Z4 


QC. Q 

-yj.o 


70 QS 

zu. yo 


1 7S 
1 . /O 


0-1-70 

y± JU 


Q-i-70 

y± JU 




1 
1 


J.L.I 


1 7-47-77 1 7 
1Z.4Z.Z / . 1Z 


,07-47- 1 7 Q 
+UZ.4J.1 j.y 


A 01 
0.U1 


70 Q 

- /u.y 


1 8 C.7 
lo.JJ 


oo 
u.uu 


1 7S-1-4S 
1 /o±40 


797-1- AO 
ZZZ± OU 


f 1 O 7Q 

nu-jy 


u 


7 1 -90 
J. l.ZU 


1 7-47-74 A 7 
1Z.4Z.J4.0Z 


_i_07-44-07 7 
+UZ.44.UZ. J 


4 71 
4. / 1 


C.4 O 
-J4.U 


1 Q 1 Q 

iy. iy 


1 70 
1. /U 


1 -i-47 
1 ±4J 


41 _i-4 A 
41 ±40 


f07 A7 
IU/-0 / 


u 


•5 1 -97 
D.Y.LI 


17-47-41 07 
1Z.4Z.41 .VJ 


_i_07-44-9C. Q 
+UZ.44.ZJ.y 


7 77 
J. / J 


77 1 
-jZ.l 


18 71 
lo.Jl 


1 80 
1.5U 


C. 1 _l7S 
-J 1±J0 


C1 _i-7S 
- J 1 ± Jo 




1 
1 


7 1 -47 
J. 1 .4Z 


17-47-71 9A 
1Z.4Z. J 1 .ZO 


, H7-4^- 1 8 7 
+UZ. 4J. 10. / 


A 1 A 
0.10 


4Q O 

-4y.u 


70 1 O 
ZU. 1U 


1 07 

i.yj 


4-1-40 
-4±4U 


4-1-40 
-4±4U 




1 
1 


7 1 -4Q 


1 7-47-77 SA 
lZ.4Z.JJ.oO 


. 07-4^-70 Q 

+uz.4j. jy.y 


C. 04 
J.y4 


47 7 
-4Z. J 


7 1 QS 

zi.yo 


1 A7 
1.0 / 


47-1-C4 
-4 / ± J4 


47-i- C.4 
-4 /± J4 




1 
1 


7 1 -An 

J.l.OU 


1 7-/17-70 77 

1z.4z.zy.zz 


_i_07-/1A-71 7 
+UZ. 40. Jl. / 


7 77 
/ . J / 


AA A 
-44.4 


1 O /17 

iy.4 / 


7 CA 
Z.jO 


Q7-1./17 

-y / ±4 / 


no 1/17 

-y /±4 / 




1 
1 


7 1 -87 


1 7-47-77 71 
1Z.4Z.JJ.Z1 


, H7-4S-10 Q 
+UZ.45. iu.y 


a 07 
o.U/ 


71 
-J l.U 


70 91 
ZU.Z1 


9 77 
L. 1 1 


70-I-C1 
-JU± J 1 


7A-1-79 
-JO± 1 L 


7 9-89 
J.Z.OZ 


1 

1 


j. 1.50 


1 7-47-78 07 
lZ.4Z.Zo.yZ 


. H7-4S- 1 A 8 
+UZ.40. 1O.0 


a 7c. 

O. / J 


7A 7 
- JO. / 


1 8 74 
lo. /4 


1 ai 

1.51 


1 7-1-78 
-LI ±J0 


1 7-1-47 
1 / ±4J 


fOQ 78 

iuy-jo 


u 


1 1 -Q7 

j.L.y I 


1 7 -A 7 -7 A 77 
1Z.4Z.jO. / / 


_l m • A Q • C 1 1 
+UZ.40.J1.1 


8 7A 
o.ZO 


97 7 
-Zj. J 


77 no 

zz.uy 


1 8/1 
1 .04 


1 A7.I.77 

-10j±jZ 


1 A7_i_79 

-10j± jZ 




i 


11-1 (Y7 


1 7-/17-77 7C 
1Z.4Z.jZ. /j 


_i_n7-/lQ-9Q A 

+uz. 4y.zy.o 


O 77 
y.LI 


97 C 
-LI .j 


1 8 8A 
1 o.oO 


n QQ 
U.oo 


1 7/1-1.97 
-1 J4±ZZ 


1 77_i-70 

-i jz±zy 


+"1 7 7Q 
11 J-Jy 


n 
U 


7 7-7 
J.Z.Z 


1 7-/17-77 1 Q 

1Z. 4Z. ZZ. 10 


_i_n7-/17-C7 7 

+UZ.4Z. jZ.Z 


7 1 n 
/. 1U 


7A Q 

- /o.y 


1 8 A1 
lo.Ol 


9 7/1 
Z. /4 


OHj./1 a 
ZU+40 


7Q-I-A 1 

jy±oi 


+" 1 n 7A 
IIU-jO 


n 


J.Z. 1 / 


1 9-/19-7C HA 
1Z.4Z.Zj.U0 


i Q7 - A 7 - /l 7 Q 

+UZ.4J.4 /.0 


A AO 

o.oy 


A7 Q 
-0 /.o 


on AQ 
ZU. Oo 


7 OQ 

j.zy 


AO-i-CQ 

-oy± jo 


AQ_i- ce 

-oy±jo 




i 
1 


1 7-71 
J.Z.Z1 


1 7-/17-7A QQ 
1Z.4Z. JO.yo 


_i_07-/l/l-07 /I 
+UZ. 44. UZ. 4 


/l 7C 
4.ZJ 


/1Q 7 

-4y. J 


71/11 
Zl .41 


7 7A 
J.JO 


9-1-/1 A 


9-1-4A 
-Z±40 




i 
1 


1 7-9A 
J.Z.ZO 


1 7-/17-77 7n 
1Z.4Z.J /.ZU 


_i_ m - /I A - 1 A 7 

+UZ. 44. 10. / 


/I 77 
4. J / 


/1A A 
-40.4 


on C7 
ZU. jZ 


7 1 Q 

j. iy 


1 O-i-CQ 

i yi jo 


1 o_i.ce 

iy±jo 




i 
1 


1 7-77 
j. L.L 1 


1 7-/17-7/1 7Q 
1Z.4Z.J4. Jo 


_i_07-/l/l - 1 Q Q 

+uz.44. iy.y 


A Q/1 
4.y4 


^1 A 
- J 1 .0 


on /ij 

ZU.4Z 


7 CC 
J. J J 


/17-uCA 
4 / ± JO 


47-i-C.A 
4 / ±J0 




i 
1 


7 7 -/in 
J.Z.4U 


1 7-/17-77 1 n 

1Z.4Z.Zj.1U 


_i_n7-/l C-Q7 Q 

+UZ. 4j.Uj.O 


7 AO 

/.oy 


An a 
-0U.4 


7 1 A7 
Zl .Oj 


7 77 
Z. 1 L 


7-1.77 
-J± 1 L 


7_i_79 
-J± 1 L 




i 
1 


1 7-/1A 

J.Z.4D 


1 7 -/1 7-7 8 On 

iz.4z.zo.yu 


_l07-/1C-91 /l 
+UZ.4J.Z1.4 


A A/1 
0.04 


C7 n 

-JZ.U 


70 QO 

zu.yu 


9 Q1 

z.y i 


1 A-x-f-A 
- 14±04 


-4±0o 


fOQ 1 7 

iuy- 1 j 


n 
u 


7 7-AA 
J.Z. DO 


1 9-49-70 97 
1Z.4-Z. JU.Z / 


4-07-4A-SQ 
-rUZ.'4-O. jy ,U 


/ . J J 


40 A 
-4U.O 


7 1 80 
Zl .oU 


9 8Q 

z.oy 


A7-t-AO 
OJ±OU 


A7-I-AO 
OJ±OU 




i 


1 7-71 
J.Z. 1 1 


19-/19-78 QC 
lZ.4Z.Zo.yj 


„l.07-/I7- 1 Q /l 

+UZ.4 / . iy.4 


a on 

o.UU 


40 Q 
-4U.O 


1 Q Q7 

iy.yz 


O 1 7 
-U. 1 J 


C-i-CQ 

-J±JO 


1 Q_uA7 

iy±o / 


+'OQ 97 

iuy-z / 


n 
U 


1 7-7C 
J.Z. / J 


17-/17-71 C7 
1Z.4Z.Z1.JJ 


_i_07-/i7-/in n 

+UZ.4 / .4U.U 


Q C/l 

y. j4 


47 Q 

-4 / .y 


79 99 
LL.LL 


O QQ 

u.yy 


1 Q-i./1/l 
- 1 o±44 


1 Q-i./1/l 
- 1 0±44 




i 
1 


1 7-27 
J.Z.oZ 


1 7-/17-77 91 
1Z.4Z.JJ.Z1 


_i_07-/l 8- 1 O Q 

+UZ.40. iu.y 


8 07 
o.U / 


71 O 
-J 1 .U 


90 91 
ZU.Z1 


9 77 
L. 1 1 


/19-i-C 1 
-4Z± J 1 


7A_u79 
-JO± 1 L 


7 1 -Q7 
J . 1 .0 / 


u 


1 7-8Q 
D.L.oy 


1 9-/19-7A 7A 
1Z.4Z.Z0. JO 


_i_07-/IQ-7/1 7 
+UZ. 40. J4. / 


Q 7Q 

y. Jo 


7Q Q 
-Jo.o 


7 1 OQ 
Zl .UO 


7/11 
J.41 


n_i_ 1 m 

U± 1 U 1 


Q 1 -i- 1 1 Q 
1± 1 10 


fOQ 41 

iuy-4i 


u 


7 7-Q7 
J.L.yL 


1 9-/19-77 /I C 
1Z.4Z. J / .4J 


_i_07-/1Q-/17 n 
+UZ.40.4J.U 


8 (\1 

o.U/ 


77 A 
-ZZ.O 


1 Q 7A 

iy. JO 


7 77 
J. JZ 


CA-i-zl 1 
-J0±41 


CA-i-zl 1 
-J0±41 




1 
i 


17-1 no 
J.Z. lUy 


1 7-/17-77 /I 7 
1Z.4Z.Z / .4J 


j_07-/1Q-78 /l 
+UZ.4y. Jo. 4 


i n ns 

1U.U0 


77 Q 
-J J.O 


70 7Q 

zu. jy 


111 
1.11 


7C-1.97 
-Zj±Zj 


S-1-40 
-o±4U 


fOQ CO f 1 7 40 

iuy-jz, iij-4u 


u 


4. 1 .J 


1 9-/19-79 QQ 

iz.4z. jz.yy 


_i_Q7 -Cn- 1/17 
+UZ. JU. 14. J 


Q Q1 

y.y i 


9C 7 
-ZJ.Z 


70 CQ 

zu. jy 


9 /17 
Z.4J 


C-1-4A 
-J±40 


C-1-4A 
-J±40 




i 

i 


4.1: 1 1 


1 O./IO-IC £fl 
1Z.4Z.jj.0U 


i no.cn. 7 1 n 
+UZ. jU.jl.y 


9.93 


-21.0 


21 .43 


0.68 


-5±53 


-J±JJ 




1 


A 1 -77 
4.1. J / 


1 9-/19-7C O/l 
lZ.4Z.Zj.y4 


_i_m-c7-77 n 
+UZ. jZ.ZZ.U 


1 7 A1 

1Z.01 


7Q 7 
-Zo. J 


7 1 77 
Zl.J / 


n A/i 
U.04 


1 Qn_i_/17 
-loU±4 / 


1 en_i./i7 
-loU±4 / 




i 


A 1 -7Q 


17-/17-7 1 /I A 
YLAL.d 1 .44 


_i_m-C7-7/i c 
+UZ. JZ.Z4.J 


1 7 nc 
IZ.Uj 


97 A 
-LLA 


ZU. / 1 


n cq 
U.JO 


1 /19-1.A7 
14Z+0Z 


1 /17_i-A9 

14Z±0Z 




i 


/I 1 -/17 
4. 1.4Z 


1 9-/19-70 Q7 

1z.4z.zy.yz 


l m-C7-/i 1 7 
+UZ.JZ.41.Z 


1 7 /l A 

1Z.40 


77 A 

-Zj.O 


1 q ec 
Io.oj 


n /i 7 

U.4 / 


1 Q9_i_97 

iyz±zj 


1 07_i-77 

iyz±zj 




i 
i 


A 1 -/1A 
4. 1.40 


1 7-/17-7Q 7 1 

i z.4z.zy . j i 


_l m - c 7 - 1 1 c 
+UZ.JJ.1 l.j 


1 7 oe 

iz.yo 


97 7 
-Z J. J 


on oo 
ZU. oo 


n OA 
u.yo 


ou± jy 


An_i-CO 

ou±jy 




i 


4.1:49 


1 O./IT-T^ CIO 

iz.4z.zj.uy 


i no ■ ci ■ o i a 
+UZ.JJ.ZJ.4 


13.61 


-27.0 


20.37 


2.94 


23±48 


23±48 




1 


A 1 -C7 
4. 1 . JZ 


1 7-/17-7/1 /1Q 
1Z.4Z. J4.4y 


_i_07-C7-7/l n 
+UZ. J J. J4.U 


1 7 QQ 

iz. oy 


17 7 
- 1 / . J 


70 A A 

zu.oo 


O 77 

u. / / 


QC-i-7/1 
0J± J4 


SC-i-74 
0j± J4 




i 

i 


4. 1 . JO 


1 7-/17-77 A/1 
1Z.4Z.ZJ. 04 


_i_Q7-C7-/1Q 1 
+UZ. J J.4y. 1 


1 A 1 A 
14. 10 


97 A 
-Z / .0 


7 1 7/1 
Zl . /4 


7 SA 
J. JO 


1 no j_/i q 

- 1UZ±40 


1 07 -i- A 8 
- lUZ±4o 




i 

i 


/I 1 -CQ 
4. I.Jo 


1 9-/19-7/1 7C 
1Z.4Z.J4. / J 


_i_Q7-C7-CQ 1 
+UZ. J J. Jo. 1 


1 7 7C 
1 J.Zj 


1 A A 
- IO.O 


1 8 8 1 
1 o.ol 


1 AS 
1 .00 


Q9_u7C 

-yz± / j 


Q7-1-7C 

-yz± / J 




i 

i 


/i 1 - An 
4.1.0U 


19-/19-78 1 Q 

iz.4z.zo. iy 


_i_n7-c/i-nc c 
+UZ. j4.Uj.j 


1 7 Q7 

i j.yz 


77 O 

-zz.y 


on A7 
ZU.O / 


7 A7 
J.O / 


1 1 7-i-C 1 
1 1 J±J 1 


1 1 7_i- C 1 
1 1 J±J 1 




i 


A 1 -AC. 
4. l.OJ 


1 9-49-9C 07 
lZ.4Z.ZJ.y / 


_i_07-C/l-9/1 9 
+UZ. J4.Z4.Z 


1 4 47 
14. 4J 


9/1 /I 
-Z4.4 


7 1 1 n 
Zl . 1U 


7 1 S 
J. 10 


77_i- /I A 
-J / ±40 


77-1-/1 A 
- J /±40 




i 
i 


A 1 -AQ 

4. i.oy 


1 7-/17-7/1 71 
1Z.4Z.Z4. / 1 


_i_07-C/l-/1/l 7 
+UZ. J4.44. / 


1 4 Q7 
14. o / 


9C. O 
-ZJ.U 


1 a AA 

lo.OO 


9/11 
Z.41 


Q7-1.C7 
0J± JZ 


S7-1-C.9 
Oj±JZ 




i 
i 


A 1 -74 
4. 1. /4 


1 9-49-78 1 A 
lZ.4Z.Zo. lO 


+UZ. JJ .U / . J 


1 A QQ 
14. oo 


91 4 
-Z 1 .4 


70 47 
ZU.4Z 


7 1 C 
J. 1 J 


1 A-1.C7 
- 10± J J 


1 A-i-C.7 
-10±JJ 




i 
i 


/I 1 -77 
4. 1 . / / 


1 7-/17-7Q 1 7 

1z.4z.zy. 1Z 


_i_07-CC-9C C 
+UZ. JJ.ZJ.J 


i c ne 

1 J.Uo 


90 1 
-ZU. 1 


1 Q SO 

iy. ju 


7 -in 

J.4U 


A1 -1-/1/1 
-Ol ±44 


A1 -1./1/1 
-0 1 ±44 




i 
i 


/I 1 -OA 
4. 1.50 


1 9 ■ /1 9 • 7 C 17 
1Z.4Z.jj.1j 


_i_m-cA-n7 q 
+UZ. jO.Uj.o 


1 C 7C 

1 j.Zj 


1 A n 
-14. U 


on /in 
ZU.4U 


9 /1C 
Z.4j 


Q9_i-CA 

-oZ± jO 


87_i- C A 
-oZ± jO 




i 


4. 1.50 


19-/19-78 CO 
1Z.4Z.Zo.JU 


_i_07 -CA- 1 A A 
+UZ. JO. io.o 


1 C Q7 

i j.yj 


1 Q A 

-iy.o 


70 81 
ZU.o 1 


1 C7 

1. JZ 


97/1-1.77 
ZZ4± JZ 


79/1-1-79 
ZZ4± JZ 




i 
i 


/I 1 -OH 

4. i.yu 


1 9-/19-77 A 7 
1Z.4Z.Z/.4J 


_i_n7-CA-77 A 

+UZ. jO.ZZ.4 


1 A 1 1 
10. 1 1 


9n 7 
-ZU.J 


7 1 1 A 

Zl. 10 


7 Q 1 
J.ol 


1 n7_i_/i 1 
-1Uj±41 


1 Q7 -1./11 

-1Uj±41 




i 
i 


4.1:97 


12:42:30.92 


+02:56:46.0 


16.21 


-17.0 


18.43 


0.90 


20±22 


20±22 




1 


5.1:4 


12:42:58.16 


+02:43:31.9 


3.07 


42.4 


21.21 


3.01 


133±26 


133±26 




1 


5.1:11 


12:42:55.57 


+02:43:53.1 


2.98 


28.5 


18.44 


1.47 


19±29 


19±29 




1 


5.1:13 


12:43:05.04 


+02:44:03.4 


4.70 


53.6 


18.88 


2.27 


-70±44 


-89±54 


f02-67 





5.1:14 


12:42:57.95 


+02:44:09.5 


3.53 


34.9 


21.71 


2.98 


70±34 


70±34 




1 


5.1:42 


12:43:02.91 


+02:45:45.0 


5.54 


36.0 


21.35 


3.05 


24±51 


24±51 




1 


5.1:51 


12:43:07.15 


+02:46:45.7 


6.99 


38.1 


20.07 


3.79 


-22±55 


-22±55 




1 



Table D.3. continued. 



ID 


R.A. 


Dec. 


^4636 


PA 


m R 


C - R 


I'spec 


Vflnal 




(J2000) 


(J2000) 


[arcmin] 


[cleg] 


[mag] 


[mag] 


[kms -1 ] 


[kms '] 


5.1:57 


12:42:58.92 


+02:47:08.1 


6.29 


21.1 


21.19 


0.87 


35±3() 


35±3() 


5.1:75 


12:43:00.47 


+02:48:15.3 


7.47 


20.7 


20.06 


2.72 


13±49 


17±93 


5.1:77 


12:43:09.11 


+02:48:28.0 


8.66 


33.7 


18.24 


-0.19 


129±21 


108±65 


5.1:81 


12:42:54.00 


+02:48:46.6 


7.58 


7.8 


22.26 


0.98 


170±52 


170±52 


5.1:88 


12:42:53.65 


+02:49:09.4 


7.95 


6.8 


19.68 


3.59 


-57±42 


-57±42 


5.1:96 


12:43:05.91 


+02:49:41.8 


9.33 


25.4 


19.82 


3.14 


-38+51 


-38+51 


5.1:99 


12:42:57.99 


+02:49:52.2 


8.84 


13.3 


20.07 


0.54 


103±52 


103±52 


5.1:104 


12:42:59.57 


+02:50:12.9 


9.27 


15.1 


21.53 


0.77 


36±65 


36±65 


6.1:2 


12:43:07.52 


+02:36:03.7 


6.82 


139.7 


21.60 


2.56 


203±43 


203±43 


6.1:5 


12:43:21.31 


+02:36:14.5 


9.32 


122.6 


19.62 


3.35 


-32+38 


-32±38 


6.1:13 


12:43:22.41 


+02:36:49.6 


9.26 


118.6 


19.49 


0.89 


-166±60 


-166±60 


6.1:18 


12:43:23.80 


+02:37:08.3 


9.43 


116.0 


21.37 


3.09 


-51 ±46 


-51+46 


6.1:22 


12:43:08.03 


+02:37:27.6 


5.92 


130.0 


19.32 


3.34 


91 ±60 


91 ±60 


6.1:35 


12:43:17.86 


+02:38:21.7 


7.57 


112.6 


21.03 


3.38 


-91±77 


-91±77 


6.1:37 


12:43:20.65 


+02:38:28.2 


8.18 


110.0 


18.34 


0.99 


36±35 


36±35 


6.1:38 


12:43:10.14 


+02:38:33.7 


5.74 


118.1 


19.45 


0.82 


-166±32 


-189+35 


6.1:42 


12:43:12.60 


+02:38:45.8 


6.20 


113.8 


20.17 


0.56 


-30±81 


-80±93 


6.1:50 


12:43:08.48 


+02:39:24.9 


5.00 


111.7 


20.21 


0.68 


-177±38 


-177+38 


6.1:53 


12:43:12.96 


+02:39:35.5 


6.00 


106.2 


20.47 


1.14 


-72±64 


-72±64 


6.1:62 


12:43:07.91 


+02:40:06.8 


4.65 


104.4 


21.00 


0.79 


-179+61 


-179+61 


6.1:75 


12:43:21.70 


+02:41:04.1 


7.95 


91.4 


21.31 


3.06 


-1()8±62 


-108±62 


6.1:86 


12:43:16.71 


+02:41:48.0 


6.72 


85.4 


21.13 


3.33 


53±72 


53±72 


6.1:87 


12:43:13.14 


+02:41:51.3 


5.84 


84.2 


20.53 


3.32 


-16±43 


-16±43 


6.1:99 


12:43:09.88 


+02:42:30.8 


5.15 


76.0 


21.52 


3.93 


-77±37 


-77±37 


6.1:100 


12:43:14.80 


+02:42:33.8 


6.36 


78.2 


22.00 


3.26 


-248±28 


-248±28 


6.2:6 


12:43:17.11 


+02:36:14.9 


8.45 


126.4 


21.89 


0.73 


3+157 


3±157 


6.2:7 


12:43:14.85 


+02:36:19.7 


7.96 


128.4 


21.81 


1.01 


-120±19 


-120+19 


6.2:9 


12:43:10.05 


+02:36:25.4 


6.99 


133.9 


18.55 


1.24 


-59±16 


-51+21 


6.2:26 


12:43:09.69 


+02:37:24.8 


6.27 


127.9 


21.81 


3.32 


44±58 


44±58 


6.2:27 


12:43:09.32 


+02:37:28.1 


6.17 


128.0 


21.88 


2.04 


-137±42 


-137±42 


6.2:50 


12:43:12.60 


+02:38:45.8 


6.20 


113.8 


20.17 


0.56 


-130±45 


-80±93 


6.2:53 


12:43:12.72 


+02:39:07.7 


6.09 


110.5 


21.29 


2.68 


23±50 


23±50 


6.2:66 


12:43:06.26 


+02:40:05.8 


4.26 


106.0 


21.43 


1.66 


-86±62 


-86±62 


6.2:67 


12:43:06.04 


+02:40:07.7 


4.20 


105.7 


18.54 


2.41 


48+81 


48±81 


6.2:77 


12:43:07.65 


+02:40:45.3 


4.47 


96.6 


19.13 


2.95 


-30±48 


-30±48 


6.2:79 


12:43:08.09 


+02:40:51.6 


4.57 


95.1 


19.82 


0.79 


-193±22 


-196±32 


6.2:84 


12:43:15.79 


+02:41:12.5 


6.47 


90.5 


20.84 


1.40 


35±32 


35+32 


6.2:93 


12:43:18.55 


+02:41:48.1 


7.18 


85.7 


21.33 


3.65 


-107±46 


-107±46 


6.2:102 


12:43:09.33 


+02:42:21.7 


4.98 


77.3 


19.24 


1.12 


89±22 


112+31 


7.1:6 


12:43:35.38 


+02:28:57.5 


16.75 


137.3 


21.44 


3.45 


-47±48 


-47±48 


7.1:7 


12:43:40.85 


+02:29:02.7 


17.65 


133.8 


20.80 


3.85 


-29±48 


-29±48 


7.1:11 


12:43:35.58 


+02:29:16.4 


16.56 


136.4 


18.65 


0.56 


-20±22 


-20±22 


7.1:15 


12:43:26.31 


+02:29:31.8 


14.85 


142.2 


19.83 


0.95 


313+20 


313±20 


7.1:19 


12:43:41.26 


+02:29:46.6 


17.23 


131.8 


21.91 


4.10 


-31 ±40 


-31±40 


7.1:20 


12:43:35.45 


+02:29:49.1 


16.14 


135.2 


20.52 


1.49 


256±27 


256±27 


7.1:30 


12:43:41.27 


+02:30:29.1 


16.76 


130.0 


20.42 


1.96 


33±34 


33±34 


7.1:31 


12:43:29.12 


+02:30:33.4 


14.52 


137.5 


21.71 


1.15 


131+35 


131+35 


7.1:33 


12:43:30.99 


+02:30:40.4 


14.75 


135.9 


21.60 


3.39 


-37±58 


-37±58 


7.1:36 


12:43:35.87 


+02:30:49.5 


15.52 


132.3 


21.92 


0.59 


70±70 


70±70 


7.1:41 


12:43:33.03 


+02:31:10.1 


14.77 


133.1 


21.03 


0.80 


260±67 


260±67 


7.1:43 


12:43:39.10 


+02:31:20.7 


15.80 


128.9 


18.54 


3.00 


-20±51 


-20±51 


7.2:5 


12:43:32.92 


+02:29:02.1 


16.28 


138.7 


22.21 


1.20 


76±199 


76±199 


7.2:19 


12:43:30.54 


+02:29:46.9 


15.33 


138.5 


21.08 


3.81 


50±85 


50±85 


7.2:31 


12:43:22.37 


+02:30:31.2 


13.47 


142.9 


22.21 


2.09 


-55±50 


-55±50 



Cross-Identif ication(s) 



t'03-45 
fi)3-47 



H2-14 
6.2:50 



f01-27 



6.1:42 



fl)2-36 



f 12-46 



Table D.3. continued. 



ID 


R.A. 


Dec. 


— p 

"4636 


PA 




C - R 


I'spec 


v final 


r^rn^^-Trlpntif ipntinn^ 1 s * 






(^J zuuu ^ 


[arcmin] 




[mag] 


[mag] 


[Km s j 


rirm c - M 

[Km s j 




7 7-Q/l 


1 7 • /I Q • If) 7 Q 
1Z.4j.3U.Z5 


+UZ. JU.4U.Z 


1 A A 7 

14. Oj 


1 ia /i 
130.4 


77 A 1 
ZZ. 01 


1 A1 

1.01 


1 0_i_A 1 

iy±oi 


1 0_i_A 1 

iy±oi 


J 


7 7-zlzl 
/ .Z.44 


1Z.4J.ZJ.UJ 


+UZ. J 1 . 1Z. 1 


1 1 ia 

1 J. JO 


1 'iS Q 
1 Jo.V 


77 1 A 
ZZ. 10 


1 'il 
1 .Jl 


ZJ±4Z 


ZJ±4Z 




/.Z.J 1 


1 0-A1-0& 70 
1Z.4J.Z0. /U 


+UZ. J 1. Jl. / 


1 J. 1 j 


1 1^ A 
1 JJ.O 


7fl Qzl 
ZU.V4 


n 7A 
u. /o 


Q7-i-An 

y / ±ou 


y / ±ou 


: 


7 7-^A 
/.Z.JO 


1Z.4J.JJ. VI 


, m.QT.nQ zl 

+uz. jz.uy.4 


1 A 78 
14. Zo 


1 70 A 

izy.o 


1 Q 7 A 
1 V. /O 


J.Uo 


R-i-zl7 
-o±4Z 


S-i- zl7 




7 7-A1 
/.Z.01 


1Z.4J.J4.ZJ 


+UZ. JZ. J 1 . 1 


1 A 1 7 
14. 1Z 


1 78 ^ 
IZo. J 


7 1 ^7 
Zl.J / 


1 aa 

1 .04 


-5l±JJ 


81 J.H 

-o 1± JJ 




7 7-A7 
/.Z.OZ 


1Z.4J.Z0. JO 


+UZ. JZ. J J.J 


1 7 Azt 
1Z. 04 


1 JJ.J 


7H nn 
zu.uu 


7 Qzl 
Z.V4 


-1 1±JJ 


-1 1±JJ 


: 


7.2:67 


12:43:33.99 


+02:32:48.7 


13.89 


127.5 


21.11 


3.28 


31 ±43 


31 ±43 




7.2:70 


12:43:28.94 


+02:32:57.7 


12.81 


130.4 


18.39 


0.00 


2±26 


2±26 




7.2:71 


12:43:32.00 


+02:33:03.0 


13.35 


128.0 


21.99 


1.76 


211+35 


211+35 




7.2:81 


12:43:32.18 


+02:33:44.1 


12.98 


125.5 


18.67 


2.89 


-21±55 


-21±55 




7.2:96 


12:43:24.95 


+02:34:33.8 


11.03 


127.4 


18.52 


0.88 


272±11 


272±11 




7.2:110 


12:43:22.13 


+02:35:22.2 


9.98 


126.2 


22.11 


2.65 


127±70 


127±70 




7.2:112 


12:43:21.07 


+02:35:27.4 


9.72 


126.7 


21.47 


0.92 


-38±22 


-38±22 





